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PEEFACE TO THE FIEST EDITION. 

Of all the celestial bodies whose motions have formed 
the subject of the investigations of astronomers, the Moon 
ha& always been regarded as that which presents the 
greatest difficulties, on account of the number of inequali- 
ties to which it is subject; but the frequent and important 
applications of the results render the Lunar Problem one 
of the highest interest, and we find that it has occupied 
the attention of the most celebrated astronomers from the 
earliest times. 

Newton's discovery of Universal Gravitation, suggested, 
it is supposed, by a rough consideration of the motions of 
the Moon, led him naturally to examine its application to 
a more severe explanation of her disturbances; and his 
Eleventh Section is the first attempt at a theoretical 
investigation of the Lunar inequalities. The results he 
obtained were found to agree very nearly with those 
determined by observation, and afforded a remarkable 
confirmation of the truth of his great principle; but 
the geometrical methods which he had adopted seem 
inadequate to so complicated a theory, and recourse has 
been had to analysis for a complete determination of the 
disturbances, and for a knowledge of the true orbit. 
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VI PREFACE. 

The following pages will, it is hoped, form a proper 
introduction to more recondite works on the subject: the 
difficulties which a person entering upon this study is most 
likely to stumble at have been dwelt upon at considerable 
length, and though different methods of investigation have 
been employed by different astronomers, the difficulties met 
with are nearly the same, and the principle of successive 
approximation is common to all. In the present work, 
the approximation is carried to the second order of small 
quantities, and this, though far from giving accurate values, 
is amply sufficient for the elucidation of the method. 

The differences in the analytical solutions arise from 
the various ways in which the position of the moon may 
be indicated by altering the system of coordinates to which 
it is referred, or again, in the same system, by choosing 
different quantities for independent variables. 

D*Alembert and Clairaut chose for coordinates the pro- 
jection of the radius vector on the plane of the ecliptic and 
the longitude of this projection. To form the differential 
equations, the true longitude was taken for independent 
variable. 

To determine the latitude, they, by analogy to Newton's 
method, employed the differential variations of the motion 
of the node and of the inclination of the orbit. 

Laplace, Damoiseau, Plana, and also Herschel and Airy, 
in their more elementary works, have found it more con- 
venient to express the variations of the latitude directly, 
by an equation of the same form as that of the radius 
vector. 

Lubbock and Pont^coulant, taking the same coordinates 
of the Moon's position, make the time the independent 
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PREFACE. VII 

variable ; and when it is desired to carry the approximation 
to a high order, this method offers the advantage of not 
requiring any reversion of series. 

Poisson proposed the method used in the planetary 
theory, that is, to determine the variation in the elements 
of the Moon's orbit, and thence to conclude the corre- 
sponding variations of the radius vector, the longitude, 
and the latitude. 

The selection of the method followed in the present-work, 
which is the same as that of Airy, Herschel, &c., was made 
on account of its simplicity; moreover, it is the method 
which has obtained sanction in this University, and it is 
hoped that it may prove of service to the student in his 
reading for the examination for Honours. In furtherance 
of this object, one of the chapters (the sixth) contains the 
physical intrepretation of the various important terms in 
the radius vector, latitude, and longitude.^ 

The seventh chapter, or appendix, contains some of the 
most interesting results in the terms of the higher orders, 
among which wiU be found the values of c and g com- 
pletely obtained to the third order. 

The last chapter is a brief historical sketch of the Lunar 
Problem before Newton, containing an account of the dis- 
coveries of the several inequalities and of the methods by 
which they were represented, those only being mentioned 
which, as the theory has since verified, were real onward 
steps. The perusal of this chapter will shew to what extent 
we are indebted to our great philosopher; at the same 
time we cannot fail being impressed with reverence for the 

* See the Beport of the << Board of Mathematical Studies" for 1850. 
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viU PREFACE. 

genius and preservance of the men who preceded him, 
and whose elaborate and multiplied hypotheses were in 
some measure necessary to the discovery of his simple 
and single law. 

I take this opportunity of acknowledging my obligations 
to several friends, whose valuable suggestions have added 
to the utility of the work. 

HUGH GODFRAY. 

8t. John's College, Cambridgej 
April, 1853. 



PEEFACE TO THE THIRD EDITION. 

lu the present edition several important alterations have 
been made, many of the parts have been entirely re-written, 
and a general re-^arrangement of the whole subject has 
taken place. In Chap. IV. Sect. 3, the necessity for the 
computation of the elements in the order ^, m, ^ is shewn. 
In the Qy>arterly Journal^ Mr. Walton pointed out the 
logical necessity for having the value of u completely correct 
to the second order before proceeding to find ^, although 
the result is the same as when the preparation for both 
values is carried on simultaneously. I have here proved 
that there is the same necessity for making the value of s 
precede that of ti, although, again, the change has no effect 
on the result. 

The term of the third order " Parallactic Inequality " will 
be found completely investigated in the appendix; and at 
the end of the volume a selection of questions from University 
and College Examination Papers has been introduced. 

Augutt, 1871. 
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LUNAR THEORY. 



CHAPTER 1. 



INTBODUOnON. 



Before proceeding to the consideration of the moon's 
motion, it will be desirable to say a few words on the law 
of attractions, and on the peculiar circumstances which enable 
us to simplify the present investigation. 

1. The law of universal gravitation, as laid down by 
Newton, is that " Every particle in the universe attracts every 
other particle^ with a force varying directly as the mass of the 
attracting particle and inversely as the square of the distance 
between themJ*^ 

The truth of this law cannot be established by abstract 
reasonings; but as it Is found that the motions of the heavenly 
bodies, calculated on the assumption of its truth, agree more 
and more closely with the observed motions the more strictly 
our calculations are performed, we have every reason to 
consider the law as an established truth, and to attribute 
any slight discrepancy between the results of calculation and 
observation to instrumental errors, to an incomplete analysis, 
or to our ignorance of the existence of some of the dis- 
turbing causes. 

Of the last source of error there is a remarkable ex- 
ample connected with the discovery of the planet Neptune, 
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2 LUNAR THEORY. 

which became known to us, as one of the bodies of our 
system,* solely by means of the perturbations it produced in 
the calculated places of the planet Uranus. These pertur- 
bations were too great to be attributed wholly to errors of 
instruments or of calculation ; and therefore, either the law 
of universal gravitation ceased to hold for a body so remote 
as Uranus, or else some unknown cause was disturbing the 
path of the planet. The first supposition was too repugnant 
to an astronomer of the nineteenth century to be entertained 
until all others had failed, — and the second supposition led 
to the detection of Neptune. The distinguished names of 
Adams and Le Verrier will be for ever connected with the 
history of this planet, and their solution of the difficult in- 
verse problem: — "Given the perturbations caused by an 
unknown planet, determine, on the assumption of the truth 
of Newton's law, the orbit and position of this disturbing 
body," — will always be considered^ one of the triumphs of 
Mathematical Science. 

Another remarkable instance of the vindication of Newton's 
law, where it seemed at first to be at fault, will be found 
' mentioned in Chap. VI. in connection with the motion of 
the Lunar Apogee. The error was owing to an incomplete 
analysis; but, whcQ the calculations were carried out more 
fully, the strongest confirmation of the law was afforded by 
the exact agreement of the results with observation. 

We shall proceed to apply the law of gravitation to in- 
vestigate the circumstances of the moon's motion ; and then 
see how it will enable us to assign her position at any time 
when sufficient data have been obtained by observation. 

2. The problem would be one of extreme, if not insur*- 
mountable difficulty, if we had to take into account simul- 

* It had been seen by Dr. Lamont at Munich, one year before its being known 
to be a planet. " Solar System, by J. R. Hind." 
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m-TRODUCtlOI^. 3 

taneonsly the actions of the earth, sun, planets, &c., on the 
moon ; but fortunately the planets are so small or so distant 
that their action may be neglected — at any rate to the order 
of our present approximation — and the attraction of the 
earth, on account of its proximity, is very much greater 
than the disturbing* action of the sun notwithstanding his 
enormous mass. We may therefore treat the question as 
that of one body, the moon, revolving about another body, 
the earth, and continually disturbed by a third body, the 
Bun. This is the celebrated problem of Three Bodies. 

Stated in this general form its exact solution has hitherto 
defied the powers of the analyst. The disturbing action 
of the sun can be expressed without difficulty in terms of 
the masses and distances, but the integration of the diffe- 
rential equations cannot be effected. If however, instead 
of retaining the whole disturbing force of the sun, which, 
as we have said, is small coptipared with the earth's action, 
we expand its expressioijxin a series and neglect all very 
small terms, it is then found possible to obtain a solution. 

This breaking up of the expression for the sun's action is 
obviously equivalent to a breaking up of the sun's force ; 
and if we afterwards wish to take into account the very 
small forces which have been thus neglected, or the small 
disturbances due to the action of the planets, the following 
principle will shew that we may do so by considering their 
effects separately, — and that the algebraical sum of all the 
disturbances so obtained will, to an order of approximation 
far beyond that which we contemplate, be the same as the 
disturbance due to their simultaneous action. We shall find 

* Since the sun attracts both the earth and moon, it is clear that its effects on 
the moon's motion relatively to the earth, or the disturbing force, will not be the 
same as the absolute force on either bodj. The absolute force of the sun on 
the moon is more than double that of the earth (Art. 99), but the relative force 
does not exceed ^^ of the earth's action (Art. 25). 
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4 LUNAR THEORY. 

an application of this principle in the investigation of the 
parallactic inequality. 

Principle of Sttperposttion of Small Disturbances, 

3. Let a particle be moving under the action of any 
number of independent forces, 
some of which are very small, 
and let A be the position of the "^ "^ 

particle at any instant. Let two of these small forces m,, m^ 
be omitted, and suppose the path of the particle under the 
action of the remaining forces to be AP in any given time. 

Let AP^ be the path which would have been described 
in the same time if m^ also had acted; AP^ differing very 
slightly from -4P, the disturbance being PP^. 

Similarly, if wi, instead of w, had acted, suppose AP^ to 
represent the disturbed path, PP^ being the disturbance, 
(^P, ^P,, AP^ are not necessarily in the same plane, nor 
even plane curves). 

Lastly, let AQh^ the actual path of the body when both 
7», and 7w, act. Join P, Q. 

The two disturbances PP^ and P^Q^ being due to the 
action of the same small force m^ on the slightly different 
paths AP and AP^^ must be very nearly equal both in 
magnitude and in direction. What difference exists between 
them must be a very small fraction of either disturbance, 
and may be neglected compared with the original path. 
We may, in fact, expect that this difference which is the 
disturbance due to m^ of a disturbance due to m, will be 
of an order compounded of the orders of the two. Thus, 
if the disturbance due to m^ and m^ be respectively of the 
second and fourth orders, the difference between P^Q and 
PP^ would probably be of the sixth order compared with 
A P. Therefore P^Q may be considered parallel and equal 
to PP.. 
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SUPERPOSITION OP SMALL DISTURBANCES. 5 

Hence the projection of the whole disturbance PQ on any 
straight line, bemg equal to the algebraical sum of the pro- 
jections of PP^ and P,Q, will be equal to the algebraical sum 
of the projections of the separate disturbances PP^^ PP^. 

Next, let there be three small disturbing forces m,, «i,, m^ 
We may consider the joint action of the two wi,, «w,, as one 
small disturbing force ; therefore, by what precedes, the total 
disturbance along any axis will be the sum of the separate 
disturbances of m^ and of the system m,, m^ ; but this last is 
the sum of the separate disturbances of m^ and m^i therefore 
the whole disturbance equals the sum of the three separate 
disturbances. 

4. This reasoning can evidently be extended to any 
number of forces ; and if or, y, z be the coordinates of the 
disturbed particle, ^(a?, y, z) any function of a?, y, «; the 
disturbance produced in (a;, y, z) will be 

u {«>, ^» ^) = ^ «* + ^ ^y + S ^^' '^^^^^s (^)*' **•» 

where Sa;=Sa;,+Sa;,+...= sum of disturbances along axis of x 

due to separate forces^ 

Sy=8y^-f 8y^+...=5 aloug axis of y, 

hz^hz^ -\-hz^+...^ along axis of is; 

therefore S0(^, y, ^) = ^S«i+^^i + ^^ 

^ dx ^^dy^^^dz ^^«^**^ 
= S,<^ (a?, y, z) 4 ijf> («, y, z) + &c., 

or total disturbance equals sum of separate disturbances, 

which establishes the principle. 

Since 4> (a;, y, z) may be the radius vector, the latitude 

or the longitude of the disturbed body, it follows that the 
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6 LUNAR THEORY. 

total disturbance in anj of these elements is the sum of the 
partial disturbances. 

5. But we must still prove another proposition, without 
which the problem, though so far simplified, would scarcely 
be less complicated than in its most general form. 

Newton's law refers to particles^ whereas the sun, earth, 
and moon are large nearly spherical bodies, and it becomes 
necessary to examine the mutual action of such bodies. 
Now, it happens that, with this law of force, the attraction 
of one sphere on another can be correctly obtained, and 
leaves the question in exactly the same state as if they were 
particles. [Princip. lib. I. prop. 75.) 

Attractions of Spherical Bodies, 

6. Let P be a particle situated at a distance OP=a from 
the centre of a sphere of uniform 
density p and radius c. The 
particle being without the sphere, 
a > c. -K-= 

Let the whole sphere be divi 
ded into circular laminae by planes 
perpendicular to OF. Let SQ be 
one of these. F8=x^ PQ = 2, and thickness of lamina = hx. 

Next, let this lamina be divided into concentric rings. 
Let R8—r be the radius of one of these rings and 8r its 
breadth, z:J?PO = 6^; 
therefore r = a? tan ^, 

Srsojsec'^^.S^. 

The attraction of an element R of this ring on the particle 

P will be -p^ along PjB, and the resolved part 

iY)QQGi OT 6l6iinpnt 
of this along PO will be 5 — 375 cos^. 
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ATTRACTIONS. 



Bnt the resultant attraction of the whole ring will clearly 
be the sum of the resolved parts along TO of the attractions^ 
of its constituent elements ; therefore, 

attraction of ring = — , ' .^ cos 5 = 27rp sin QZx . 8^ ; 
•17 sec (7 

therefore, attraction of whole lamina SQ^^trpZxx sin O.dO 

Jo 
=27rpSa;fl--V 

Again, «• = oj* + c* — (a — a;)' = 2ax — (a* — c*) ; 
therefore zSz — aSx^ 

and attraction of lamina = 2'irp ( —^ — 8zj ; 

/. attraction of whole sphere = 27r/} \- — | - ^-w-2— \ 

(from z=:a — c to z = a + c) 

_ 47rpc' 
"■ 3a' 

where Jf = mass of sphere = — ^ . 

o 

Hence, the attraction of the whole sphere is precisely the 
same as if the whole mass were condensed into its centre. 

Cor. 1. The attraction of a shell (radios c and thickness 
he) will be obtained from the preceding expression by diflTe- 
rentiating it with respect to c, and is 

,^ ^. p , „ 47rpc'. he mass of shell 
attraction of shell = — ^ — = i , 

the same as if the mass were collected at its centre. 
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8 LUNAB THEORY. 

Cob. 2, Therefore, the attraction of a heterogeneous sphere 
on an external particle will be the same as if the whole mass 
were condensed into its centre^ provided the density he the same 
at all points equally distant from the centre^ for then the 
whole sphere may be considered as the aggregate of an 
infinite number of uniform shells. 

7. Let us now consider the case of one sphere attract- 
ing another. Suppose P in the preceding article to be an 
Elementary particle of a sphere M\ whose centre 0* suppose 
at a distance a from 0. Then, since action and reaction 
are equal and opposite, P will attract the whole sphere M 
just as it would do a particle of mass M placed at 0. The 
same is true of all the elementary particles which compose 
the sphere M\ therefore the sphere M* will attract the sphere 
J/ as if the whole mass of the latter were condensed into its 
centre ; but the attraction of the sphere M* on a particle 
is the same as if the attracting sphere were condensed into 
its centre 0' ; therefore, 

Two spheres attract one another as if the whole matter of 
each sphere were collected at its centre. 

8. This remarkable result, which, as may be shewn, holds 
only when the law of attraction is that of the inverse square 
of the distance, or that of the direct distance, or a com- 
bination of these by addition or subtraction, reduces the 
problem of the sun, earth, and moon to that of three par- 
ticles. The slight error due to the bodies not being perfect 
spheres is here neglected, being of an order higher than 
that to which we intend to carry the present investiga- 
tion : this error, however, though very small, is appreciable, 
and when a nearer approximation is required, it becomes 
necessary to have regard to this circumstance. (See Ap- 
pendix, Art. 109). 
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CHAPTER 11. 



MOTION RELATIVE TO THE EARTH. 

9. When a number of particles are in motion under their 
mutual attractions or other forces, and the motion relatively 
to one of them is required, we must bring that one to rest 
and then keep it at rest, without altering the relative motions 
of the others with respect to it. 

Now, first, the chosen particle will be brought to rest 
by giving it at any instant a velocity equal and opposite to 
that which it has at that instant ; secondly, it will be kept 
at rest by applying to it accelerating forces equal and 
opposite to those which act upon it. 

Therefore give the same velocity and apply the same 
accelerating forces to all the bodies of the system, and their 
absolute motions about the chosen body, which is now at rest, 
will be the same as their relative motions previously. 

Problem of Two Bodies. 

10. As the sun disturbs the moon's motion with respect 
to the earth, it is important to know what that motion would 
have been if no disturbance had existed, or generally : — 

Two todies attracting one another with forces varying 
directly as the mass and inversely as the square of the dis- 
tance^ to determine the orbit of one relatively to the other. 

Let Mj M! be the masses of the bodies, M being the body 
whose motion relatively to if is required, r the distance 
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10 LUNAR THEOEY. 

between them at any tune t^ and the angle between r and 

some fixed prime radios. 

M 
The accelerating force of if on M' equals -7 acting towards 

M^ while that of Jf' on JIf equals —^ in the opposite direction. 

Therefore, by the principle above stated, we must apply to 
both M and M* accelerating forces equal and opposite to this 
latter force, and M* will move about Jf fixed, the accelerating 

force on M* being — -^ — = )tAtt*, if fi = M-\-M' and r = -. 
Hence, -jg^ +11 = ^^, 

where A = r* -^ = twice the area described in a unit of time, 

and Integrating 

u = ^,{l+6Cos(^-a)}, 

e and a being constants to be determined by the circumstances 
of the motion at any given time. 

This is the equation to a conic section referred to its 

focus, the eccentricity being e, the semi-latus rectum — , and 

the angle made by the apse line with the prime radius a. 

By observations made on the actual path of the moon, 
which cannot differ very widely from its undisturbed orbit, 
we infer that this latter would be an ellipse with an 
eccentricity about ^. 

In the same manner, if the sun and earth be the two 
bodies considered, the relative orbit would be an ellipse with 
an eccentricity about -^q. 

11. The angle 5- a between the radius vector and the 
apse line is called the true anomaly. 
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Two BODIES. IX 

If n is the angular velocity of a radius vector which, 
moving uniformly, would accomplish its revolution in the 
same time as the true one, both passing through the apse 
at the same instant; then nt + B-a is called the mean 
anomaly where e is a constant depending on the instant 
from. which the time is reckoned, its value being the angle 
between the prime radius and the uniformly revolving one 
when ^ = 0. 

Thus, let MT be the fixed line or prime radius, 

A the apse, 

M* the moving body 
at time t^ 

MfA the uniformly re- 
volving radius at same time, the 
direction of motion being repre- 
sented by the arrow. 

And, let MD be the position of Jf/i when f = 0, 

then tMD = e and is called the epoch,* 

DMfA^nt^ 

tMA = a = longitude of the apse ; 
therefore, mean anomaly = AMii, = w^ -f e - a, 
trueanomaly = uliOf' = YiOf'-TJf^ = «-a. 

12. To express the mean anomaly in terms of the true in 
a series ascending according to thepoioers of Q^as far as e*. 

27r ^2 area 

= 2^- 




periodic time A 

2irh h 

'2^^ a^^Jil'-ey 



* We shall see (Art. 86) that the introdnctioa of the epoch e is avoided 
in the Lunar Theory bj a particular assumption ; but in the Planetary it forma 
one of the important elements of the orbit. 
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12 LUMAB THEOBT' 

therefore A=no*(l — ^)*, 

dt r* a*ll-e*f 1 



M h h * {l.+eco8(^-a)j* 

= ^(l-e')*{l + ecos(^-a)P 

= - (1 - f e") {1 - 28 co» (^ - o) + 3^ COS* (^ - a)}: 

= - {l-2eco8(^-a) + f«»co82(^-a)}; 
n 

therefore n^ + e = 5 — 2« sin (^ — a) + |e' sin 2 (^ — a), 

or (w« + e-a) = (^-»)-26sitt(^-a) + ife"8m2(^-a), 

the required relation. 

13. To eocpreaa ike true anomaly in terms of the mean to* 
the same order of approximation. 

tf-a = we+e-a + 2esin(^-a)-^"8in2(^-a)e..(l); 

.•. ^ — a=s/ie + 8-a first approximation. 

Substituting this in the first small terms of (I), we get . 

5 — a = w^ + e — a + 26 sin (n^ + 8 — a).-a second approximation* 

Substitute the second approximation in that small term^ 
of (]) which is multiplied by e, and the first approximation 
in that multiplied by e*, the result will be correct to that 
term, and gives 

5 - a = w^ + 8 — a + 26 sin {n^ + 8 - a + 2e sin {nt-\-e — a)} 

- f e* sin 2 [nt + e - a) 

= n< + 8 - a + 26 wn (wf + 8 — a) 

+ 46* cos(n<+s-a) sin (w«+8-a)— 16' sin2(7iHe- a) 

= n< + 8 - a+ 2e sin(n< + s - a) + f e* sin2 («^ + 8 - a), 

the required relation. 
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THREE BODIES. 13 

The development could be carried on by the same process 
to any power of e, the coefficient of ^ would be found 

■J^fsinS (w<+8— a) -J sin (w« + e — a), 

but in what follows we shall not require anything beyond e*. 

Prollem of Three Bodies. 

14. In order to fix the position of the moon with respeet 
to the centre of the earth, which, by means of the process 
described in Art. (9), is supposed brought to rest, we must 
have some determinate invariable plane passing through th& 
earth's centre to which the motion may be referred. 

If the sun and earth were the only bodies in the universe^ 
then the plane of the elliptic orbit (in which, according 
to the last section, the motion would take place) would be 
fixed, and might be taken as the plane of reference; but,, 
as soon as we take into account the action of the moon 
and planets — especially the moon — the plane ceases to be 
fixed, and some other plane must be found not affected by 
these disturbances. 

Theory teaches us that such a plane exists,*^ but as its 
exact determination can only be the work of time, the 



♦ See Poinsot, ^^Theorie a determination de Vequateur du systeme solaire^ 
where he proves that an invariable plane exists for the solar system, that is, a 
plane whose position relatively to the fixed stars will always be the same^ 
whatever changes the orbits of the planets may experience; but as its p($sition 
depends on the moments of inertia of the snn, planets, and satellites, and there- 
lore on their internal conformation, it cannot be determined a prion, and ages 
must elapse before observation can famish sufficient data tot doing so a 
posteriori. 

This result Poinsot obtains on the supposition that the solar system is a 
free system ; but it is possible, as he furthermore remarks, nay probable, that th& 
stars exert some action upon it, it follows that this invariable plane may itself be 
variable, though the change must, according to our ideas of time and space^. 
be indefinitely slow and smaU. 
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14 LUNAR THEOEr. 

following theorem will supply us with a plane whose motion 
is extremely slow, and which may for a very long period^ 
and to a degree of approximation far beyond that to which 
we shall carry our investigations, be considered as fixed and 
coinciding with its position at present 

In what follows the action of the planets b too small 
to be taken into account. 

15. The centre of gravity of the earth and moon describea 
relatively to the sun an orbit very nearly in one plane and 
elliptic; the square of the ratio of the distances of the moon 
and sun from the earth being neglected.* 




Let 8j JSj M be the centres of the sun, earth, and moon, 
O the centre of gravity of the last two. Now the motion 
of O is the same as if the whole mass ^+ if were collected 
there and acted on by forces equal and parallel to the moving 
forces which act on E and M. The whole force oil O \b 
therefore in the plane 8EM\ join 8Q. 

Let L 8GM^ ft), and let m' be the sun's absolute force. 
Moving force on E^ ■ '^ in E8j 

moving force on M^ '^ in M8. 



* This ratio is about 7^, and, as we ehall see Art. (21), such a qnantitj 
we shall consider as of the 2nd order of small quantities, and its square therefore 
of the 4th order. Our investigations are carried to the 2nd order only. 
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PLANE OF REFERENCE. 15 

These applied to O parallel to themselves are equi- 
valent to 

— fi^ 'SW — ^ direction GM^ by the 

I triangle 

m\E.8G m'.M.SG ^^. ^f forces, 

ana ^^^ + ^r^s "^aj, 



but 



E M^ _ M+E 
GM^ GE'^ ME ' 



whence E. GE^ M. GM^ (2f + E) ^^^ • 

Therefore substituting, and dividing by if+-S,2^e get 
accelerating forces* 

, GM.GE ( 1 1 \ . J- r mur 

Again, 8E*^8G'+ GE'-j-^SG.GEcosto 



= 8G' ( 



, 2(yJS; . GE*\ 

l + __eosa, + -g^.j; 

ZGE 



therefore -^^ = -^^ ( 1 — ^75- cosw )| 

SE"" bG \ bG J\ omitting quantities 

. ., , 1 1 /, ZGM \\ to be neglected. 



* In Btrictnefis it would be necessary, sinoe we have brought S to rest to • 
apply to both M and E^ and therefore to (?, accelerating forces equal and opposite 
to those which E and M themselves exert on 8 ; but the mass of /9 is so large 
compared with those of E and M^ that we may safely neglect these forces in 
this approximate determination of the path of G^ the error being of a still 

~sg) ' 
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16 LUNAR THEOBY. 

Therefore the acceleratiDg force in the direction OM 
Zm'.GMME 



8G' 



cosa> 



= - 3 (accelerating force of sun on G) -^^ . -^^ cosco 

= 0, according to standard of approximation adopted. 
And the accelerating force in G8 

m\8G {E-^M B(E.GE-M.GM) ) 

^MVEts^ m ^'^} 

t 

= "o^oi t^ ^^ same approximation. 

Hence, the force on G^ is a central force tending to S 
and varying inversely as the square of the distance ; there- 
fore the orbit of G about S is very approximately an ellipse 
with 8 in the focus; and the plane of this ellipse is, as far 
as our investigations are concerned, a fixed plane. 

This fixed plane is called the plane of the ecliptic^ or 
fiimply the ecliptic. 

16. We shall now suppose a plane through the earth's 
centre parallel to the ecliptic; and this, when the earth is 
brought to rest, will give us the fixed plane of reference 
we require (14); the ecliptic then making small monthly 
oscillations from one side to the other of our fixed plane. 

17. Since G describes an ellipse relatively to the sun, 
the sun will describe the same ellipse relatively to G ; but, 
as seen from the earth, the orbit of the sun will be slightly 
different, as we shall now shew. 

First : The sun will have a latitude, that is, will be out 
of the plane of reference ; and this latitude will be of the 
same name (north or south) as that of the moon, and 
deducible from it when the ratio of the distances of the two 
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ERRORS IN sun's POSITION. 17 

bodies from the earth, and that the masses of the earth 




and moon are known. For if S'EAT be the fixed plane 
through E^ and S\ 0\ if' the projections of 5, (?, Jf, then 
SES' is the sun's latitude, and MEM' is the moon's; and 
the two bodies are obviously both on the same side of 
the plane. 

w 7ii.x S8' GG' EG. ain(moon's\^t) 
sm(sun's lat.) = ^^ - ^g^ = -^ 

M EM . , , , , 

Now, it is known that M is about ^ of jS, 

and^Jf :jJo oiE8] 

- . - , 1 x\ sin (moon's lat.) , 

therefore smfsun s lat.) = — ,^,^, nearly. 

^ d240U •' 

And as the moon's latitude never exceeds 5*" 9', the sun's 
latitude will always be less than 1". 

Secondly, as to the sun's longitude: let -E^ be the 
direction of the first point of Aries, — ^that is, a fixed line 
in our plane of reference from which the longitudes of the 
bodies are reckoned. Then TES* = ff is the sun's longitude 
as seen from J?; but, as seen from G^ the sun's longitude 
would be the angle between G8 and a line through G (not 
drawn in the figure) parallel to Ef. 

The difference between these two longitudes is the angle 
E8'G'] tind, it EG' 8' = a>, 



sinE8' G* =^^ sm (o = 32405 approximately ; 
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18 LUNAR THEORY. , 

therefore sinES'G' never exceeds ^^i^^, and E&O' is a 
small angle not exceeding 7". 

Also E8' ^S'0'< EG' < ^^ij^jj S' G\ 

18. Now, by assuming the longitude and distance of the 
sun as seen from E to be the same as when seen from G^ 
we commit the above small errors in the position of S] 
that is, we assume the sun to be at S" instead of 8^ 8' 8" 
being drawn equal and parallel to G'E. If our object were 
the determination of the sun's position, it would be neces- 
sary to take this into account ; but when investigating the 
motion of the moon, the disturbing action of the sun will 
not, on account of his great distance, be appreciably altered 
by supposing him to be at 8" instead of 8; and the errors 
thus introduced are much too small to be considered. 

Hence we may assume that the motion of the sun about 
the earth at rest is an ellipse having the earth for its focus, 
and its equation 

w' = a'{li-6'cos(e'-?)}, 

where f is the longitude of the sun's perigee and e' = ^ is 
the eccentricity, and we are safe that no appreciable error 
will ensue in the determination of the moon's place.* 

♦ That is, as far as the three bodies alone are concerned ;— but, since the 
attractions of the planets may, and in fact do, disturb the elliptic orbit of 
the sun about (7, the same cause will disturb the assigned orbit about E, A 
remarkable result of this disturbance is noticed in Appendix, Art. (108). 
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CHAPTER III. 



RIGOROUS DII*FERENTIAL EQUATIONS OP T?HE MOON*S MOTION 
AND APPROXIMATE EXPRESSIONS OF THE FORCES. 

19. The earth having been reduced to rest by the process 

described In Art. (9), let us take its ^ 

centre E for origin. Let r, be the IZ^^^^"^ 

polar coordinates of the projection M' "^^^^^^^^^^^^ 

of the moon on the fixed plane of re- ^ 

ference, d being the longitude TEM' measured from a fixed 

line j^ in that plane. Let s be the tangent of the moon's 

latitude MEM\ so that MM' = ra. 

Next, let the accelerating forces which act upon the moon 
be resolved into these three : 

P parallel to the projected radius M'E and towards the earth, 
T parallel to the fixed plane, perpendicular to P, and in the 

direction of 6 increasing^ 
8 perpendicular to the fixed plane and towards it. 

The equations of motion will be 

;a(40=^---v B- 

S-'(f)'=-^- -P). 

-de ^ ("'^- 

20. These three equations for determining the moon's 
motion take the time t for independent variable ; but it will 
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20 LUNAR THEORY^ 

be more convenient in the following process to consider the 
longitude as such, and our next step will be to change the 
independent variable from t to 0. 
From (i) we get 

dt dt\ dt)~ dt' 
therefore (r' ^= A' + 2 [Tr^dO, 

= IF suppose, whence E -jg= Tr\ . 

Therefore ^^ = -« = Evf. if w = - , 

at r ^ r ' 

an important eqttatton connecting, the time t with the hngi-^ 
tude 6. 

, . dr dr dO „du 

atr at \ au/ cUf \ <w/ 

Substitute in (ii), 

therefore Hu' ^ (E^^ + S'u' = P, 



or 



^»'&"+»)+»-K-^- <^'f 



and substituting for E' and for E^ their values h'+ ^hd0' 

T . . . 

.and -5 , then dividing by A'u', and transposing, 

this is called the differential equation of the moon^s radius 
vector^ 
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DIFFERENTIAL EQUATIONS. 21 



Lastly, 



a ^\u) dl ^'dtf Ue de\_d (jjf da du\\ 
TYt i( d^8 d*u\ TjdH , / ds du\ 

from (A P. = fiv(,g4-t.)+^fu4;; 
therefore i%- 5= jff'u' (^ + A+^^u* J; 



tborefore 

cTs Ps-S T d» (d^8 \ [T ^^ 

dff*'^''''l^^''¥^'d»'^W*'^Vm'' 



this is the differential equation of the moorCa latitude. 

21. K these three equations eould be integrated under 
these general forms, the problem of the moon's motion would 
be completely solved ; for as only four variables u, 0^ a^ and 
t are involved (the accelerating forces P, T^ and 8 being 
functions of these), the values of three of them, as «i, d, a^ 
eould be obtained in terms of the fourth t\ that is, the radiu» 
vector, longitude, and latitude would be known correspond- 
ing to a given time. 

But the integration has never y«t been effected, except 
for particular values of P, T, and 8 ;< and the method which 
we are in consequence forced to adopt is that of successive 
approximation, by which the values of Uy 0j and a are 
obtained in series in which the terms proceed according to^ 
ascending powers of certain small fractions. Some one- 
fraction being chosen as a standard with which all other 
are compared, the order of the approximation is esteemed 
by the highest power of the small fractions retained. 
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22 LUNAR THEORY. 

It is usual to consider ^^ ^^ ^ small fraction of the first order, 

consequently ^ of ^V^^iir ^^ second 

^T^TF^s third 

and so on, other fractions being considered as of the 1st, 
2nd, &c. orders, according as they more nearly coincide with 
iki ^hj5i &c. 

22. It is necessary therefore, before we can approximate 
at all, that we should have a* previous knowledge (a rough 
one is sufficient) of the values of some of the quantities 
involved in our investigations; and for this knowledge we 
must have recourse to observation. 

We shall therefore assume as data the following results 
of observation : 

(1) The moon moves in longitude about thirteen times 
as fast as the sun. The ratio of the mean motions in longi-* 
tude represented by m is therefore about 3^^, and may be 
considered as of the 1st order.* 

(2) The sun's distance from the earth is about 400 times 
as great as the moon's distance. 

Hence the ratio of the mean distances =* ^f^ is of the 
second order.f 

(3) The eccentricity e^ of the elliptic orbit which the 
sun approximately describes about the earth is about ^, 
and this, approaching nearer in value to ^ than to ^^o, 
will be considered as of the 1st order. 

(4) During one revolution, the moon moves pretty ac- 
curately in a plane inclined to the plane of the ecliptic at 

* This approximate value of m is easily obtained ;~thA moon is foond to 
perform the tour of the heavens, returning to the same position among the 
fixed stars, in about 27| days; the sun takes 865^ days to accomplish the 
same journey. 

fThe distances of the luminaries may be calculated from their horizontal 
parallaxes, found by observations made at remote geographical stations. (See 
the Author's Astroiwmy^ chap. xvi.). 
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EXPRESSIONS OP THE PORCES. 23 

an angle whose tangent is about 3^, and therefore of the 
Ist order.* 

(5) Its orbit in this plane is very nearly an ellipse having 
the centre of the earth in its focus, and whose eccentricity is 
about equal to our standard of small fractions of the 1st order, 
viz. ^; and this will also be very nearly true of the projec- 
tion of the orbit on the plane of the eclipticf 

To calculate the values of the forces. 

23. We are now in possession of the data requisite for 
beginning our approximations, and we shall proceed to the 
determination of the values of P, T, and S in terms of the 
coordinates of the positions of the sun and moon. 

Let 8j Ej M be the centres of the sun, earth, and moon. 



m\ E^ M their masses, 

E\ M\ the projections on the plane of the ecliptic, 

O the centre of gravity of E and M. 

* That the moon's orbit during one revolution is very nearly a plane inclined 
as we have stated, will be found by noting her position day after day among the 
fixed stars; and the sun's path having previously been ascertained in a similar 
way, the rules of Spherical Trigonometry will easily enable us to verify these 
statements. {Astronomy ^ p. 92). » 

t The elliptic nature and the value of the eccentricity of the moon's orbit may 
be found by daily observation of her parallax, whence her distance from the earth's 
centre may be determined : corresponding observations of her place in the heavens 
being taken, and corrected for parallax to reduce them to the earth's centre, will 
determine her angular motion. Lines proportional to the distances being then 
drawn from a point in the proper directions, the extremities mark out the form 
of the moon's orbit. 

A similar method applied to observations of the diameter of the sun will 
determine the eccentricity of its orbit. {^Astronomy, p. 138). 
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LUNAR THEORY. 



The forces we have to take into account are, the forces 
which act directly on Jtf, and forces equal and opposite to 
those which act on E; — these last being applied to the whole 
system so that £J may be a fixed point. 



Forces on M, 



Attraction of E^ 



E 



.mME 



.in ME 



ME'" 
attraction of 8] m!.MG 

. 1 .X I cm^B par* to gig 
^qmvalent to ^ SM '^ 



Forces on E. 
M 



Attraction of Jlf= 



attraction of 8 
equivalent to 



MB'" 
m'.EG 



An EM 



.in EM 



BE* 
-^j^ pai^.to QB 



Therefore, the whole attraction upon Jf, when E is 
brought to rest, is 

E+M^ ,fMG ^EG\. „_, 

and m'SG (ot?? *" *o^) parallel to GS. 

These are the rigorous expressions of the accelerating 
forces on M^ and they can be expressed in terms of the 
masses and coordinates of the bodies. 

Approximate values of P, 37, S, 
24. Our investigations will be carried only to the second 
order ; it will be sufficient therefore if, in the preceding, we 
neglect small quantities of the fourth and higher orders. 
Let/i* = ^+2f, 

MGM' = tan"*5 = moon's latitude, 

^'T the direction of the first point of Aries, 

56^ = r' = A ; ^rE'8= ff = longitude of sun, 

M'E = r = i ; iTCE'M = tf = longitude of moon, 
/. iSB'if' = ^—^ = difference of longitude of sun and moon. 
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EXPRESSIONS OP THE FORCES. 25 

Now, SJir = jgitf'»+20f'» 

= 80^+ OM'* - 28G. OM' cos SQM' + MM* 



r"(l- 



^GW „^.., GM\ 
2 — r- COS SGM +—^] ; 



therefore -^, = /» {^ + ~y~ cos (^ - ^j| ; 

for e-ff or fi'j;'Jf' diflfers from 8GM' by less than 7', 

Art. (17), aod ( — j-) is neglected, being of the fourth order. 

Similarly, -g^, ~ "^ ] * ~ cos(^ - ^)l ; therefore 

the accelerating forces on the moon are approximately 

-^ + ^ {MG + GB) in direction ME, 

and ^{GM'+ GE') co8(5- &) paraUel to G8] 

T 

whence P= {-^ + ^' M^ cosMGM'- ^ ME" co8'(5-^) 

= ^^^ r+-5^--5-5 {l+C0S2^-^)} 

= /i«''(l-i«»)-^{i + |cos2(d-5')l, 
T= - ^' M-E' cos (5 - ^) sin [6 - ff\ 

= -|2^%In2(5-n 

S=(^ + ^ilf2?)8iniffl'Jf 

- f M otVV(14»') ] « . 



P, _ 5= _ ^^ {| + 1 cos2 (^ - *■)}. 
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26 LUNAR THEORY. 

The differential equations in Art. (20), when these values 
of the forces are substituted in them, would contain the 
variables vi ^d ff\ but v! is given in terms of & by the 
equation of Art. (18) w =a'fl +«' cos(^ — f)}, and we shall 
find means to establish a connexion between Q and ff^ which 
will enable us to eliminate the latter. 

25. Before proceeding further, it will be important to 
consider the order of the disturbing effect of the sun*s action, 
compared with the direct action of the earth. If we examine 
the values of P, T, and 8^ it will be found that the most 
important of the terms containing w', which are clearly the 
disturbing forces since they depepd upon the sun, are in- 
volved in the form —^ , while those independent of the sun's 

action enter in the form -s . 

We must therefore find the order of — ^ compared with -^ , 

orof-7i -|. 

r r 

Now the orbits being nearly circular, and m the ratio 
of the mean motions. Art. (22), we have 

_ mean motion of sun periodic time of moon 
~~ mean motion of moon periodic time of sun 

therefore -^:^::m* : 1, 

T T 

or the disturbing force of the sun is of the second order 
when compared with the direct action of the earth. 
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CHAPTER IV. 

INTBGBATION OF THE DIFFERENTIAL EQUATIONS. 

Section I. 

General process described. 

26. If in the differential equations we substitute the 
values of P, 1\ and 8 just obtained, we are led to the 
following forms: 

which cannot be integrated ; but it is found possible to obtain 
a solution which will be a first approximation, if we retain 
only the important terms on the right-hand side of each 
equation. This first approximation being substituted in the 
equations, the terms next in order of importance may be 
retained and a new solution becomes possible, which in its 
turn is the basis of a further approximation, and so on. 

First, neglect the disturbing force of the sun which is 
of the second order, and also the moon's latitude, which, 
as will be seen by referring to the expressions for the 
forces (Art. 24), will either enter to the second power or else 
in combination with the disturbing force. 
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28 LUNAB THEOKY. 

When this is done the equations become ihtegrable, and 
values of u and 8 may be obtained in terms of 6 correct to 
the first order; this valae of u will then enable us to get 
the connexion between 6 and t to the same order. 

[Let OS, however, bear in mind that the equations- thus integrated are not 
the differential equations of the moon's motions, but only approximate forms^ 
of them; and it is, therefore, possible that the results obtained may be ap- 
proximations to the true solutions only for a limited range of valuer. 

Whether they are so or not, will appear by comparing them with what we- 
already know of the motion from obsenration; and this previous knowledge^ 
in the event of their not being approximations, will probably suggest such, 
modifications of them as will render them so.] 

To proceed to the next approximation :' — We must sim- 
plify the right^iand members of our general equations by 
employing the values of u and a just obtained, retaining: 
now terms of the second order which we bad previously 
neglected. 

The integration can always be performed when the right- 
hand members of the equations are circular functions (sines 
or cosines) of or of its multiples; and as our first approxi- 
mation will give us the values of u and s in that form, the 
new expressions for the forces will also be of the same 
character, and the equations be again integrable. 

Thus we shall obtain new values of u, «, t correct to the* 
second order. These values, introduced in the same manner 
in the second members and terms of the next higher order 
retained, will lead to a third approximation, and so on, to^ 
any order; except that if we wish to carry it on beyond 
the third, the approximate values of the forces, giv^n in 
Art. (24), would no longer be sufficiently exact, and we 
must obtain more correct values from the expressions o£ 
Art. (23). 

27. We must here stop to notice a peculiarity in these 
equations, when solved by this process. We have said that 
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GENERAL PKOCESS OF INTEGRATION. 29 

to obtain the values to any order^ all terms up to that order 
must be retained in the second members ; but it may happen 
that a term of an order beyond that to which ^e are work- 
ing would, if retained, be so altered by the integration as 
to come within the proposed order. 

Such terms must therefore not be rejected| and we 
shall proceed to examine by what means they may be 
recognised. 

28. Suppose then that after an approximation to a certain 
order, the substitutions for the next steps haye brought the 
equation in u to the form 

^ + 11= -K?cos(j>^+^)+ , 

where the coefficient Q is one order beyond that which we 
intend to retain. The solution of this equation will be of 
the form 

tt= + ff' cos(/?5 + a')+ 

O* being a constant which may be determined by substi-^ 
luting this value of « in the differential equation, 

whence Q' = i ; 

from which we learn that if p differs very little from 1, O' 
will be at least one order lower than Q^ and will come 
within our proposed approximation; and consequently the^ 
term O co&(p0-¥ H) must be retained in the differential 
equation. 

The equation of the moon's latitude being of the same 
form a» that of the radius vector, the same remarks apply 
to it. 

29. If p is small, 0'=0 very nearly, and the tern> 
G COB {p6 + H) will not rise in importance in the value of u. 
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30 tUNAR TliEORY» 

We should therefore be led to reject it in the diflferentlal 
equation, except for the following reason : — 

In finding the connexion between the longitude and the 
time (one of the principal objects of the Theory), we must 
use the equation 

Now, substituting for u its value in terms of 0y we shall 
find that such a term as O' cos(/?^-f JT) in u gives rise to 
a term O" cos(p0 + H) in this equation, where O' and O' 
are of the same order, then 

J= G"coi,{p0-\-H), 

whence <=« — Bin{pd'\-H); 

therefore, when p is a small quantity of the first order, — 

will be one order lower than G'\ and the term will have 

risen in importance by the integration. 

T 
But yet further, if such terms occur in n-s j they will be 

twice increased in value ; for they increase once in forming 

Ittz ^^1 A^d ^^^^ again, as above, in finding t, 

30. We have, therefore, the following rule : — 

In approximating to any given order to the values of u 
and Sj we rnust^ in their differential equations^ retain periodical 
terms ONE ORDER beyond the proposed one^ when the coefficient 
of 6 in their argument^ is near unity. 

♦The angle of a periodical term is its only variable part and is called 
the argument. 
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TERMS TO BE RETAINED. 31 

In approximating to the valite of f, we must^ in the same 
equations^ also retain terms ONE ORDER beyond the proposed 
one when the coefficient of the argument is near zero ; and 

T 

when such terms occur in p— , they must he retained TWO 

ORDERS heyond the proposed approximation.^ 

31. We shall here, for convenience of reference, bring 
together the equations and the approximate expressions for 
the forces : — 

d^u P T du ^(d^u \ fT ,^ , ., 

'8 T ds ^fd's \ t T j^ ,T^. 



d^ _Ps^8 
dd''^^" h'u' 

T.-wMi-^'T («). 

^. =f.(I-M- ^ll + 3c<«2(9-«')) iD), 

^ = -IW-^I*-") <^' 



* Instead of the forces which really act on the moon, we originally substituted 
three equivalent ones, P, T^ 8) these again are, by the preceding expressions, 

p 
replaced by a set of others. For, we may conceive each of the terms in tj— j , Ac, 

to correspond to a force, — a component of P, T, or 8; each force having the 
same argument as the term to which it corresponds, and therefore going through 
its cycle of values in the same time. Now, by Art., (29), when the coefficient of 
in the argument is near unity, the term becomes important in the radius vector, 
and when near zero, in the longitude : hence, a force whose period is nearly the 
same as that of the moon, produces important effects in the radius vector ; and a 
force whose period is very long will be important in its effect on the longitude. 
See Airy's Tracts, Planetary Theory, p. 78. 
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32 LUNAR THEORY. 

Section II. 
To solve the Equations to the first order, 

32. Neglect all terms depending on the disturbing forcei 
t.6. those which contain m*\ such terms being of the second 
order. Art. (25). 

The latitude s of the moon can never exceed the inclina- 
tion of the orbit to the ecliptic ; but this inclination is of the 
first order, Art. (22), therefore s is at least of the first order 
and I? may be neglected. 

Whence -, = a*; j;i = 0; —^ =0, 

and the differential equations become 



cTu /jk 



de* 



+ u=n 



d^s 

whence w = p {1 + e cos (fl — a)} ; or, writing a for ^ , 

tt = a {l + ecos(^-a)} (CTj), 

and « = Ajsin(^-7) v (/SJ, 

«, a, A;, 7 bemg the four constants introduced by- integration. 

33. These results are in perfect agreement with what 
rough observations had already taught us concemmg the 
moon's motion Art. (22); for tt = a{l+6C0s(^-a)} repre- 
sents motion in an ellipse about the earth as focus. 

And 5 = A; sin (d— 7) indicates motion in a plane inclined 
to the ecliptic at an angle tan~*Ar. 

For, if TOM' be the ecliptic, 
M the moon's place, MM' an arc 
perpendicular to the ecliptic, then 
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SOLUTION TO FIRST ORDER. 33 

^M'^0] and if TO be taken equal to 7, and OM joined. 
by an arc of gteat circle, we have 

tan J/Jlf ' = tanif Oil/ ' sin OJIf ' ; 

of 5 = tSLnMOM' sin (5 - 7), 

which, compared with the equation above, shews thslt 

MOM' = tiin% 

Therefore the moon moves in a plane passing through 
a fixed point and making a constant angle with the 
ecliptic. 

The Values of a and 7 introduced in the above solutions 
are respectively the longtitude of the apse and of the node. 

34. What the equations cannot teach us, however, is 
the magnitude of the quantities e and k. For this we must 
have recourse to observation, and by referring to Art. (22), 
we see that e is about ^ and k about -j^, that is, both 
quantities are of the first order. Their exact values cannot 
yet be obtained t the means of doing so from multiplied 
observations will be indicated further on. 

35. Lastly, to find the connexion between t and 0, the 
equation (C) becomes, r=0, 

d^ J^ ^ J 1 

de ^ hu* " ha' [l + e coa{0 - a)\' * 

!Now this is the very same equation that we had found 
connecting t and in the problem of two bodies^ Art. (12), as 
we ought to expect, since we have neglected the sun's action. 
Therefore, if p be the moon's mean angular velocity, we 
should, following the same process as in the article referred 
to, arrive at the result 

^=^^ + e + 26sin(^ + e-a) + fe*sin2(;?<-f e-a)H- , 

which is correct only \o the first order^ since we have 
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34 ^.UNAR THEORY. 

rejected some terms of the second order by neglecting the 
disturbing force. 

36. The arbitrary constant e, introduced in the process 
of integration, can be got rid of by a proper assumption: 
this assumption is, that the time t is reckoned from the 
instant when the mean value of is «ero.* 

For, since the mean value of fl, found by rejecting the 
periodical terms, is pt + e] if, when this vanishes, < = 0, we 
must have e = ; therefore 

^pt-\-2e 8in(j[>t -- a) 0^, 

correct to the first order.f 

37. We have now obtained three results, ITj, 8^^ 0^, 
as solutions to the first order of our differential equations, 
and we must employ them to obtain the next approximate 
solutions; but, before they can be so employed, they will 
require to be slightly modified — in such a manner, however, 
as not to Interfere with their degree of approximation. 

The necessity for such a modification will appear from 
the following considerations : 

* When a function of a variable oontains periodical terms which go through 
all their changes positive and negative as the variable increases continuously, 
.the mean value of the function is the part which is independent of the periodical 
terms. 

t We shall also employ this method of correcting the integral in our next 
approximation to the value of in terms of t ; and if we purposed to cany our 
approximations to a higher order than the second, we should still adopt the 
same value, that is, zero, for the arbitrary constant introduced by the integration. 
To shew the advantage of thus correcting with respect to mean values : suppose 
we reckoned the time ^m some definite value of 0, for. instance when 6 = 0; 
then, in the first approximation, 

= « + 2« sin (e - o) 
is the equation for determining the constant c, and in the second approximation, 
s would be found from 

= « + 2e sin (e - o) + |e« sin 2 (« -o) + , 

giving different values of • at each successive approximation. 
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MODIFICATION OF FIRST SOLUTION. 35 

Suppose we proceed with the values already obtained; 
we have, by Art. (24), 



,'8 



J 



= a+ H--4 cos(fi- «)+. 

and this being substituted in the differential equation (^A) 
of Art (31;, gives 

^ + M = a+ + u4cos(^-a)+ , 

the solution of which is 

w = a {1 + e cosC^ - a)} + + \AB sin(^ - a). 

Our first approximate value w = a {1 + e cos(5 — a)} is here 
corrected by a term which, on account of the factor ^, admits 
of indefinite increase, and thus becomes ultimately more 
important than that with which we started. Such a cor- 
rection is inadmissible, for our first value is confirmed by 
observation (22) to be very nearly the true one. The moon's 
distance, as determined by her parallax, is never much less 
than 60 times the earth's radius; whereas this new value 
of u^ when 6 is very great, would make the distance in- 
definitely small. On the same principle, we see that any 
solution, which comprises a term of the form -4^ sin (^ — a), 
cannot be an approximate solution except for a small range 
of values of 6. 

Such terms, * if they really had an existence in our system, must end in its 
* destruction, or at least in the total subversion of its present state; but when 
*they do occur, they have their origin, not in the nature of the differential 
' equations, but in the imperfection of our analysis, and in the inadequate repre- 
^sentation of the perturbations, and are to be got rid of, or rather included in 
< more general expressions of a periodical nature, by a more refined investigation 
<than that which led us to them. The nature of this difficulty will be easily 
'understood from the following reasoning. Suppose that a term, such as 
^a sin (^0 + ^), should exist in the value of «<, in which A being extremely 
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36 J.UNAR THEOKY. 

* minute, the period of the inet^oallty denoted by it wonld be of great length ; 
' then, whatever might be the value of the coefficient a, the inequality would 

* still be always confined within certain limits, and after many ages would returti 

< to its former state. 

'Suppose now that our peculiar mode of arriving at the value of u led us 

* to this term, not in its real analytical form a sin {AB + B)y but by the way of 
*fts development in powers of 0,'o + /30 + ]^6' + Ac. ; and that, not at once, but 

* piecemeal, as it were; a first approximation giving us only the term a, a 

* second adding t;h§ term /30, and so on. If we stopped here, it is obvious that 
'we should mistake the nature of this inequality, and that a really periodical 

< function, from the efEect of an imperfect approximation, would appear under 

'the form of one not periodical These terms ip. the value of u, when they 

'occur, are not super^uous; they are essential to its expression, but they lead 
'us to erroQeops ponciusions f» to the stability of our system and the general 
'laws of its perturbatioiis, unless we l^eep in view that they are only parts of 

* series; the principal parts, it is true, when we confine ourselves to intervals 
^ of moderate l^ngth, but which cease to be so aftier th^ lapse of very long timet 
' the rest of the series acquiring ultimately the preponderance, and compensating 
' the want pf periodicity of its ^rst terms.-SiR Jo^N Hb^scihbl, Encyclopadia 
Metropolitctnqr-TsjBlOAL ABf RONOMY, p. 679. 

38. To e^tritcate ourselves from this difficulty, and 
to alter tbe solution so that none but periodical terms 
may be introduced, let us again observe tbat the equation 

-j^ + v=p = a, which gave the solution Z7, and thus led 

to the difficulty, is only aA approximate form of the first 
order of the exact equation. Any value of m, therefore, 

which satisfies the approximate equation -;j^ + u=^a to the 

first order, and which evad,es the difficulty mentioned above, 
may be takeu as a solution to the same order of the exact 
equation. 

Th^se conditions will be satisfied if we assume 

w = a {1 + 6 cos(c^ - a)} ; 

then ^, + w = a + 06 ( 1 - c') cos(c^ - a) 

= a to the^first order, 
provided I — c' be of the first order at least. 
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MODIFICATION OP FIBST SOLUTION. 37 

39. The introduction of the factor c in the argument is 
an artifice due to Clairaut; and its effect, as we shall see 
in Art. (66) is equivalent to supposing that the apse of the 
moon's orbit is not fixed. Now, although the observations 
recorded in Art. (22) did not suggest this, we must bear in 
mind that they were extremely rough, and carried on only 
for a short interval; but when they are made with a little 
more accuracy, and extended over several revolutions of the 
moon, it is soon found that both her apse and the plane 
of her orbit are in constant motion. 

Clairaut was fully aware of these motions, and there is 
no doubt that he was led to the above form of the value 
of u by that consideration, and by his acquaintance with the 
results of Newton's ninth section, which, when translated into 
analytical language, lead at once to the same form.* 

We might, therefore, taking for granted the results of 
observation, have commenced our approximation at this step, 
and have at once written down u = a[l-\-e copied — a)} ; but 
we should, in so doing, have merely postponed the difficulty 
to the next step, since there again, as we shall find^ the 
differential equation is/of the form 

— ^ + u = a function of fi, 



* Newton l^as there shewn, that il the angular Telocity of the orbit be to 

G^ — F* 

that of the body as G^ - F to (?, the additional centripetal force is - — -^ — **'**» 

the original force being ijlu\ Therefore 





■4- 


(P 


«, 






F* fi 


F^ 



F 

where 7; is the same as our e. 



u = ah +€008^— 6-oj|, 
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38 LUNAR THEORY. 

the correct integral of which would be 

u = Acos(0-'JS) + , 

and this would at the next operation bring in a term with 
6 for a coefficient, which we now know must not be. We' 
shall, therefore, hereafter omit such terms as A coh{0^B) 
altogether, and merely write 

u = a{l + €C03(c^- a)} + 

40. The value of c will be found more and more cor- 
rectly at each successive approximation, by always assuming 
these terms as the first two terms of the value of u, then 
substituting in the differential equation and equating coeffi- 
cients. It will thus be found that ae (1 — c^) must equal the 
coefficient of cos(cd — a) in the differential equation ; and this 
will enable us to determine c to the same degree of approxi- 
mation as that of the differential equation itself. See 
Arts. (52) and (103). So far, all that we know about c is that 
it differs from unity at most by a quantity of the first order. 

41. In carrying on the solution of 5, the same difficulty 
arises as in u^ and it will be found necessary to change it into 

s = k sin( gO -y)j 
g being a quantity which differs from unity at most by a 
quantity of the first order. 

The introduction of g is connected with the motion of the 
node in the same way as that of c is with the motion of the 
apse (85) ; and the value of g will be determined by a process 
similar to that explained above for c. Arts. (51) and (102). 

42. The connexion between 6 and t will also be modified 
by this change in the value of w, 

dt_^ 1 1 

dd ^ ha^ {1 + 6 cos(c6^ - a)j" * 
d.ct 1 1 



or 



d.cd Aa* {1 + e cos(c^ - a)}' ' 
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MODIFICATION OF FIRST SOLUTION. 39 

Here cd and ct hold the places which and t occupied 
in (35); therefore 

c0 = cpt+ 2e sin (cp* - a), 
or =pt + 2e 8in(cp« — a), 

to the first order, since - = e to the first order. 

43. Since the disturbing forces are to be taken into 
account in the next approximation, we shall have to use 
the value of v! found in (18), which is 

t£' = a'{l+6'cos(^'-?)}- 
but this introduces ff ; we must therefore further modify it 
by substituting for 0' its value in terms of fl, and it will 
'be found suflScient, for the purpose of the present work, 
to obtain the connexion between them to the first order, 
which may be done as follows : 

Let ra be the ratio of the mean motions of the sun and 
moon, 
p\ p their mean angular velocities ; .\p = mp^ 

pt-\-l3jpt mean longitudes at time «, 13 being the sun's 

longitude when * = 0, 

0\ true longitudes at time f, 

f, a longitude of perigees when ^ = ; 

therefore 0' — 5*= sun's true anomaly, 

and p't + 13''^^ mean anomaly. 

But, by Art. (13), 
true anomaly = mean anomaly + 2e' sin (mean anomaly) + &c. ; 

therefore ^' = pt-¥ + 2e' sin {p't + ^ - f ) + 

= mpt-^/S +2e'&m[mpt+0-^) + 

= m0 + +2e'sin(m^ + i8-f) 
to the first order ; 
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40 LUNAR THEORY. 

because jpt^^d — ^e sin (cd - a) to the first order by (42). 
Whence w' = a {1 + e' cos (tnd + /8 — f)} to the first order. 

44. The values of sin 2 (^- ^) and cos 2 (0-0') can also 
be readily obtained to the same order : 

sin2(fl-^) 

= sin {(2 - 2m) ^ - 2;8 - Ae' sin {m0-\-fi- ?)} 

= sin{(2-2m)^-2;8}-4e'sin(w^+;8-g')cos{(2-2m)^-2/8} 

= sin {(2 -2m)0^ 2)8} -2e' sin {(2 - m) ^ - )8 - ?} 

+ 2e' sin {(2 - 3m) - 3^8 + f }. 

Similarly, co82 (^ - fi') = cos {(2 - 2m) ^ - 2)8} 

-2e' co8{(2 -m) ^- )8- ?} + 2e' cos{(2 - 3m) ^- 3/3+ ?}. 

The first term of each of these is all we shall require. 

Section III. 

To solve the equations to the Second Order. 

45« Let us recapitulate the results of the last approxi-^ 
mation. 

tt = a {1 + e cos(c^ - a)}, 
tt' = a' {1 + e' cos(mfl + )8 - 1;% 
8^k8in(g0-'y\ 

0^ff^{\^m)0- ^'•2e'«Ln{m0-¥fi-Ky 

These values must now be substituted in the expressions 
for the forces, retaining terms above the second order, when, 
according to the criterion of Art. (30), they promise to be- 
come of the second order after integrating. 
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SOLUTION TO THE SECOND ORDER. 41 

The differential equations [A) and {B) will then assume 
the forms 

and the integration of these will enable us to obtain u and s 
to the second order; after which, equation '{C) 'will give the 
conuexlon between 6 and t to the same order. 

46. The quantity -ytt , which we shall meet with as a 

coefficient of the terms due to the disturbing force, can be 
replaced by m'a, m being the ratio of the mean motions of 
the sun and moon. 

So long as we neglected the disturbing force, h and a 
had determinate values : — they belonged to the ellipse which 
formed our first imperfect solution, and would therefore be 
known from the circumstances of motion in that ellipse at 
any instant; h being double the area described in a unit 
of time, and a the reciprocal of the semi-latus rectum. It 
would consequently be impossible to assume any arbitrary 
connexion between them. But, when we proceed to a second 
approximation and introduce the disturbing force, there is 
no longer a determinate ellipse to which the h and a apply : 
the equation /i = A*a of Art. (32) merely shews that a and h 
must refer to some one of the Instantaneous ellipses which 
the moon could describe about the earth if the disturbance 
were to cease, and we are at liberty to select any one of 
these which will allow us to proceed with our approximation. 
The particular ellipse is determined by the above assumed 

relation -vira =^i*«> and the selection is suggested and 

O 
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42 LUNAR THEORY. 

justified by the following reasoning : 

_-P' _ average period of moon about earth 
p periad of sun about earth ' 

but, since the instantaneous ellipses are nearly circles, we 
have, as in Art. (25), 

(period of moon about earth)' _ w'a" , 
(period of sun about earth)" "" /ita' ^ ' 

therefore if a be properly chosen, 



o 772 a 



'" "" AV • 



47. If we examine the equations of Art. (31) it will be 
seen that we cannot solve the differential equation in u to 
the second order until we know the value of « accurately to 
that order ; for, suppose 

« = A;sinC^5-7)+ 7,, 

where V^ contains the terms of the second order not yet 
found; and suppose that among the terms in V^ there be 
one i/sin(Z^ — X), then this value of s substituted in the 
expression {B) gives 

^ = j^.{l-f[Asin(;y^-7) + isin(Z^-X)...7}-... 

= a {1 - P' sin^C^fl - 7) - ZUL sin [gO - 7) sin(?^- X)...}, 
and equation {A) becomes, so far as this term is concerned, 

+ |A:i[cos{(<7 + 0^-(7 + X)}].... 

Now, if it should happen that Z is a small quantity, then 
the two coeflScients g-l and ^ + Z will be both near unity, 
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SOLUTION TO THE SECOND ORDER. 4S 

and therefore these two terms, although of the third order 
in the differential equation, would become important in the 
value of u. Again, if I is near 2, then l-g will be near 
unity, and the same remark will apply. 

Such terms do not really occur in the value of «, but 
there is none the less a necessity for verifying the fact. 

If I is near unity, then the coefficient g — l becomes small ; 
and although the corresponding term in u is still of the third 
order, it becomes of the second order, and therefore impor- 
tant in e. Art. (30). 

If I is not near 0, 1 or 2, the term will not affect, to 
that order, the values pf u or of t. 

48. Again, it is necessary to have the value of u correct 

T ' 
to the second order, before we can compute rr^ to the fourth 

order, with a full assurance that we have secured all terms 
of that order in which the coefficient of the argument 
is nearly zero*; — such terms becoming, as we have seen 
Art. (30), important in the value of t 
For, assume that 

tt = a{H-6Cos(c^-a)+ TTJ, 

^, T 37nV»{l + 6'cosfw^+)8-?)}' . ,.. _ /^ _^. 

then Tr-% = - nTTTTi 7-71 ^ r^ sin 1(2 - 2w) 5- 2^} 

A*m" 2AV {1 + e cos [06 - a) + fF,|* ^^ ' ^^ 

=- f m" (1...- 4 TTJ sin {(2 - 2m) 5 - 2^9}. 

Therefore if W^ contain a term R cos [r6 — p) of the second 
order where r is near 2 (except 2-2wi) there will be a 



* This was pointed out by Mr. Walton in the Quarterly Journal of Mathematics, 
Vol. iz., p. 227; and in this edition the text had been modified in accordance 
with his suggestion ; but the method given by him for finding the value of u is 
defective on account of its not mentioning the necessity for the previous deter- 
mination of ». 
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44 LUNAR THEOKT. 

T 

corresponding term in p— 3 of the fourth order with a small 

coefficient of 0^ and therefore an important term in the 
value of U 

49. For the sake of simplifying the expressions we shall 
generally omit a, ^8, and 7 in the arguments. If we remark 
that c and a always appear together in the form cO — a, it 
will be sufficient to write c5, instead of cd — a; and, for 
the same reason, we may replace md-\-l3 and gff — y by 
mO^ and gd^. The suffix (J will remind us of the omission, 
and we may at any stage re-introduce the omitted symbols. 



To compute s to the second order, 

50. Since we must first determine the value of 5, it 
becomes necessary to examine the differential equation (JB) 
to see whether any difficulty of a similar kind to those we 
have been here considering may arise, or whether the 
elements on which the computation depends are already 
known with sufficient accuracy. We fortunately find that 
all the terms of the right-hand member of the equation are 
of the third order, and therefore that the values already 
obtained will ensure accuracy to the second order, and we 
may safely begin our work here. 

51. To the first order we have 

da . /IT A (because ^ = 1 

^"*^'^''^^'-*«»''^''' {to first order, 



^+-0. 
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TO COMPUTE 8 TO THE SECOND ORDER. 45 

Ps^ 8 
In the expression for , 3 , which is of the third order, 

it will be sufficient to put a! for w , and a for m, then 
"AV*"^ f-^^Asin^^,{l + cos(2-2m)5,} 

= - f Tw'A sin^r^j + fw'^A sin (2 - 2w -5^) ^i, 

the other term of the third order being rejected because the 
coefficient of is not near unity. 

The next term of equation {B) is rr-a-jni ^°^ ^^ In '^^ 

T 
of the first order, it will be sufficient to have t5— « correct 

to the second order, 

^» = -Ksin(2-2m)^., 

T ds 
■j^ . ^ = - \m^k sin (2 - 2m) 0^ cos^^^ 

= - \m*k sin (2 - 2w - g) 0^^ 

neglecting the other term of the third order because the 
coefficient of is not near unity. 

r T 
which =0 to the third order, since Itt^^^ depending on 

the disturbing force is of the second order, and the other 
factor = to the first order. 



The equation in s becomes 

— |m'A sin^^j 

+ |?n'A;8in{2-2w-^)^j. 



To solve it, assume 



hAng0^ 
+ ^sin(2-2m-flr)5^. 
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46 LUNAR THEOET. 

Sobstitate, and equate coefficients of like terms, 

h {\-g') = -lm^k', .: g^{l + f m'ji = 1 + f m', 

Therefore, accurately, to the second order,* 

a-h mL[g6 - 7) 4 Imk sin {(2 - 2?w - jy) ^ - 2^8 + 7}.../8',. 

The term of the second order, having a coefficient near 
unity, need not be taken into account in finding the value 
of u (47) ; but it will have to be brought in subsequently 
to determine t. 

To determine u to the second order. 

62. We must compute the right-hand member of the 
equation (-4), Art. (31), accurately to the second order, and 
include those terms of the third order which have a coefficient 
near unity. 

To the first order we have 

w = a(l + e cosc6'j), 

du . ^ . ^ (because c = 1 

ad * » (to first order, 

From the expressions (2>) and [E)^ Art. (31), we get 

^, = a(I-fA»sinVJ 

mV {J + e' cos(m^, - f )}' . „ ,^ « n /i > 
- 2AV(l + ecoscV^ {1 4 3 C08(2 - 2«) ^J, 

♦ No complementary term P cos(0- Q) is added j for, though by the theory 
of differential equations, this would form a necessary part of the solution,, we 
haye seen Art. (39) that it cannot, in this shape, form a port of the correct yalne^ 
but will be comprised in the terms whose argument is ^0- y. 



Digitized by 



Google 



TO DETERMINE U TO THE SECOND ORDER. 



47 



- \m* {1 + 3e' cos [m0^ - ?) - 3e cosc^,} {1+3 cos (2 - 2m) ^J, 

1 P 

- ^rir = 1 - 1*' + P' cos 25r^j - f m* - |m* cos (2 - 2m) ^, 
ci /» i£ 

+ |m*e cosc^j + |m'e cos (2 — 2m - c) ^j. 
The other terms of the third order, which arise in the 
development of the expression, are neglected because the 
coefficients of are not near unity, 

T 3m'a"{l + e'co8(m^,-?)}» . . ^ . . 

= - f m' [1 + 3e' cos(m^, - C) - 4e cosc^J sin (2 - 2m) 6^ 
= - |m^ sin (2 - 2m) 6^ + 3m*e sin (2 - 2m - c) 6^^ 

^, 8 • 3^ = f m*ae sin (2 — 2m) ^^ sinc^^ 

= f mV cos (2 — 2m — c) ^, ; 
rejecting the other terms of the third order according to (30). 

TT-s ^^=1 cT-TT- COS (2 - 2m) 6^ - - — ^ — cos (2 - 2m - 6)0^ 
h^u"" ^2-2m ^ ^ * 2-2m-c ^ ^ * 

= |m' cos(2 - 2m) 0^ - 3m*e cos(2 - 2m - c) ^j, ' 

2 (^ +M^ Lr-, rf5=|m'a cos (2-2m) 0,-6m'ae cos(2-2m-c)^j. 

The equation in u becomes 

+ f m*e cosc^j 

4- |A' cos 2^5, 

- 3m' cos (2- 2m) ^, 

+ ym'e cos(2 - 2m - c) 0^. 



h 



To integrate this, assume 

u 

a 
l+A 
+ ecosc5, 
•f J? cos 2^^, 
+ Ocos(2-.2m)^, 
+ Z)cos(2-2m-c)^i. 
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Substitute, and equate the coeflScients of like terms : 



C{1-(2-2»i/}=-3ot', 



15m*e , - 



i>= 



2{1-(1-27W+|7WY} 



w = a 



.K 



i> {1 - (2 - 2m - c)'} = ^/mV 
Therefore 

+ e cos (c^ — a) 
-iA:' cos 2 (^6^-7) 
+ 7w*cos{(2-2m)5-2/9} 
, + V^«cos{(2-2m-c)^-2;8 + a} 
This is the accurate value of \i to the second orden 

To find t to the second order. 
53. The equation to be considered is 

and, making use of the values of a and u obtained to the 

second order, we must retrace our steps to recompute u and 

T 
YT-s so as to include those terms of the third order in u in 
n u 

which the coefficient of the argument is near zero, and 

T 
similar terms to the fourth order in p-. Art. (48). To do 

this we need not repeat the whole work, but only so much 
of it as will enable us to pick out the additional terms in 
question. 
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TO FIND t TO THE SECOND ORDER. 49 

-im*{l+3e' co8(m^ -0-36 cosc5,l{l+3 C08(2-2wi)dJ 
= l.-.-fwA' siDj^^j sin (2- 2m- ^) 5^- f wiV coH(m6^-0 
= l...-^mA' cos(2 - 2m - 2g) 0^ - fmV C08(m^j- (;), 

r 3m^a^ (1 + e' cosmg^ - ?y sin (2 - 2m) g, 

AV"" 2A'a* (1 + e coscV, - JA* co82^^J* 

= - f m* sin (2 - 2m) ^, {1 + e co&c0^ - i^j' cos2^^J^ 
= - f m« sin (2 - 2m) ^, {(1 + e cosc^,)"^ + k' cos2g0^} 
= - f m* 8in(2 - 2m) e, j.-lOe* cos^ctf^ + A:' cos25r5J 
= - f m* sin (2 - 2m) 0^ {...5e* cos2c0j^ + k' eoB2g0^] 
= - V^«* 8in(2-2m-2c)^j-fmT 8in(2-2m-2flr) 0^. 
These are the only important additional terms. 

j^d0=^- ^vni cos (2- 2m -2c?) 0-\inV cos(2-2m-2flr) ^^, 

because 2-2m-2c=-2m+fm', and 2-2m-2^=s-2m-fm\ 
Substittrting in the differential equation for u, the new 

rp J 

terms in 75-5 ^^ will be of the fourth order and must be 

-j^ + »j \if-^d0 tnH bring in the two 

terms 

- ^me^a cos (2 -2m-2c) 0^ and -fmi'a cos(2 -2m-2gr)^j 5 
and we get, so far as important additional terms are eoncemed^ 

- fe + ttl^ 1...+ V^e* cos(2 -.2m-2c) 0^ 

+ (- X + f)»^<«>s(2-2m-25r)^-|wVcos(m^j-?), 
and the corresponding terms in u will also be of the third 
order, viz. 

- =... ^me" cos(2 - 2m - 2c) 0^ + -^mk* cos(2 - 2m - 2g) 0^ 

- f mV cos (mtfj — ?). 
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50 LUNAR THEORY. 

The complete value of u therefore, so far as importaut terms 
are concerned for determining the value of f, is 

- = 1 - |A' - ^m' + e cosc^j - ik' cos2g6^ 

+ w' cos(2 - 2m) 0^ + ^me cos (2 - 2w - c) 0^ 

+ *^oi* cos (2 - 2m - 2c) 0^ + ^mh^ cos (2 - 2«n - 2g) 0^ 

— §7wV cos(m5j — f), 

and the complete value of \^jj-^ d0 will be 

\ji-^d0=^\m* coa(2-2m)0^-^^m' cos(2- 2m-2c)^, 

- |«wAj* cos (2 - 2w - 25^) 5j, 

the other term of the third order has no influence on tj the 
coefficient of not being small. 

54. We are now in a position to find the value of « to 
the second order : 

If- = l + ^coso5, + TF where TF contains all terms of the 
second and third orders in the expression above, then 

= ^, {I - 2 (€ coso^^H- TF) + 3 (e cosc^,+ TF)'- ...} 
= ^, (1 - 2e cosc^j - 2 W+ 3e" cos*c5, + GeTFcosc^,)... 
= ^ (1 + |e' - 2e cosc^,- 2 W+ f e* cos2c^,+ 6e TFcosctf^). 

y Google 
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TO FIND t TO THE SECOND ORDER* 61 

The only part of GeTTcosc^, which gives a term of the third 
order in which the coeflScient of the argument is small, is 

6ecosc^j {...+ ^me cos (2- 2m — c)^^ ...}, 
of which one term is 

+ ^me^ cos(2 - ^m - 2c) 6^. 

Therefore ^i is equal to 

(1 + |e' + %k^ + m'-'2e coscd^ + f e* co82c0^ 
+ Ik' cos2g0^ - 2m* cos(2 - 2m) 0^ 
- ^me cos (2 - 2m - c) ^, + 3m V cos(m^j - f) 
- i/me* cos (2 - 2m - 2c) 0^- f m^' cos (2 - 2m - 2g) 0^. 
Also from above 

rr ri--|m*cos (2-2m) 0^+^me' cos (2-2m-2c) 0^ 

J Vu"" { + \mW cos (2 - 2m - 2g) 0^. 

We have now to multiply these results together, and we see 

that the terms having for coefficients ^m^ and |m^* will 

disappear in the product.* These were the terms which 

T 
were originally of the fourth order in y^-, . 

fl + f e' + f A' + m^ - 2e cosc^i + f e* cos2c^i 
H- \k^ cos2g0^ - V^* cos(2 - 2m) 0, 
- ^mc cos(2 - 2m - c) ^j + 3mV cos(m^j - ?). 



♦ The principle of the superposition of small disturbances might have led as 
a priori to anticipate the disappearance of these terms: — for, when approximating 
to the second order we may, according to this principle, calculate separately the 
effect of the disturbing action of the sun, which is of the second order, and that 
of any other ind^enderU disturbance of the sam^ or of a higher order, and then 
add the results. No term therefore which bears a trace of the action of both 
these forces should present itself to the second order. 

Now since neither k nor e enters into the argument of any term, it is obvious 
that a term whose coefficient contains k^ or c* at any stage will retain them 
through aU subsequent operations ; and if the terms having {2 — 2m — 2g) 6, 
or (2 — 2w» — 2c) 6i for argument could remain in the result to the second order 
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52 LUNAB THEOBT. 

Let ^,(i + |,« + K + m')=i, 

therefore 2> = Aa' (I — f e* — |A* - vr?) to the third order ; 
therefore, multiplying by p and integrating, we get, still to 
the second order, 
jp«=5-2e sin(ctf-a)+|e* sin2 (cO-a) 

- V»w' sin {(2 - 2wi) ^ - 2^8} 

- V'wesin{(2-2m-c)5-2i8+al 
+ 3rWsin(?n^ + ^-f) 

DO constant is added, the time being reckoned from the 
instant when the mean value of vanishes, for the reasons 
explained in (Art. 36). 

Coordinates expressed in terms of the time. 

65. The preceding equations D^, /S,, 0^ give the reci- 
procal of the radius vector, the latitude and the time in 
terms of the true longitude ; but the principal object of the 
analytical investigations of the Lunar Theory being the 
formation of tables which give the coordinates of the moon 
at stated times, we must express u^ 5, and 6 in terms of t. 

To do this, we must reverse the series pt==0 — &c., and 
then substitute the value of in the expressions for u and s. 

Now ^ =sp< + 2£ sin (c0 - a) to the first order 

=^pt-\-2e sin(cp^ — a) ; 

fhsy Tfoold th^r@ have k^ or e? r^p^ctively 'in the coefficient, since they haye 
them wheQ they ^t appear. 3nt B, e^ and m being independent, such terms 
most haye arisen from the combination of others which originally (in the 
expressions of the forces) involyed these quantities separately: — terms like 
4;^ cos 2 (^0 — y) and e* cos 2 {c0 — q) combined with terms inyolving m : that 
is, terms which are representatiyes of disturbing forces of the second order (see 
liote p. 81) combined with others depending on the sun's action. Such terms 
must therefore be of an order beyond the second. 
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therefore c^ — a = cp« — a + 26 sin (cp< — a) to the first order, 

2e Bin(cff — a) = 2e {sin (cp< —a) +26 sm(cp^— a) cos(cp^— a)} 

to the second order) 

= 26sin(c5p« — a) + 26* sin 2 (cp^- a) ; 

and as and pt differ by a quantity of the first order, they 
may be used indiscriminately in terms of the second order ; 
therefere 

^=2?^f 26 sin (cpt — a)+ f 6* sin 2 (cpt — a) 

— JA* sin 2 (ffpt — 7) 

+ ^m* sin {(2 - 2m)pt-2l3} ... ©'., 

+ ^me Bm[(2'-2m-c)pt-2^+a] 

— Sme' sin(7np^ + ^ — 

56. In the value of u given in Art. (52), substitute pt for 
6 in terms of the second order, and pt + 26 sin(cp< - a) in the 
term of the first order ; then 

1— |A;*-im*— 6*+6cos(cp<-a)+6* cos2(cp<-a) "^ 
"ik* coH2(gpt- 7) 
+ m' cos {(2 - 2m) pt - 2^8} 
+ ^me cos{(2-2m-c)^«-2i8+a} 



U=iW 



.u;. 



57. Similarly, the expression for a becomes 
^ = i sin [(jjpt — 7) + 26 8in(cp^ - a)} 

+ fm* sin {(2 - 2m - ,y) jp« - 2i8 + 7} ; 
'Ar\(()pt-i) 
+ 6 8in{(^ + c)p<-a-7} 
— 6 sin \(g — c)pt + a — 7} 
+ §w sin {(2 - 2m - ^);?f ^ 2i8 + 7} 

The expression for 8 is more complex in this form than when 
given in terms of the true longitude 0. 



or 



8^k 



s;. 
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54 LUNAR THEORY. 

MoorCs Parallax. 

58. If P be the moon's mean parallax, and U the parallax 
at the time f, 

' radius of earth R r» ^- t .n . ,. , , 

-v/(l+0 

= -Btt {1 — JA' 4- iJ(^ cos2 ((^pf -^ 7)} to the second order, 

II - k' — ^m' - e* + 6 coa(cpt — a)+e* cos2 (cpt- a) 
+ wi* cos {(2 - 2m) ;?/ - 2/8} 
+ i^5wiecos{(2-2m-c)^«-2/8 + a} ; 
but P= the portion which is independent of periodical terms, 

il + e cos(cp« - a) + e' C082 (cpt — a) 
+ 7n* cos {(2 - 2m) pt - 2/3} 
+ V^6 cos {(2 - 2m - c) j?^ - 2/9 + a} 
neglecting terms of the third order. We see that, to the 
second order, the variable part of the parallax is independent 
of the inclination* 

59. Here we terminate our approximations to the values 
of w, 5, and 0. If we wished to carry them to the third 
order, it would be necessary to include some terms of the 
fourth and fifth orders according to Art .(30), and the ap- 
proximate values of P, T^ and 8^ given in Art. (24), would 
no longer be suflSciently accurate, but we should have to 
recur to the exact values, and from them obtain terms of 
an order beyond those already employed.* The process 
followed in the preceding pages is a sufficient clue to what 
would have to be done for a higher approximation. 

* See Parallactic Inequalitt/y Art. (105). 
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SOLUTIONS TO A HIGHER ORDER. 56 

The coordinates u and & of the sun's position are, by 
the theory of elliptic motion, known in terms of the time f, 
and t is given in terms of the longitude d by the equation 
0,. Hence vi and & can be obtained in terms of d\ but 
it will be necessary to take into account the slow progressive 
motion of the sun's perigee, which we have hitherto neglected. 
This we may do by writing dff — f for ff — f, c' being a 
quantity which diflfers very little from unity.* 

These values of u\ d\ together with those of u and a 
given by TJ^ and fi^„ are then to be substituted in the 
corrected values of the forces, and thence in the differential 
equations. The integrations being performed as before will 
give the values of w, a, and t in terms of 6 to the third 
order, and from these, as in Arts. (55), (56), and (57), may be 
obtained u^ a, and 6 in terms of t. 

60. More approximate values of c and g are obtained at 
the same time, by means of the coefficients of cos {cQ — a) 
and sin {gd — 7) in the differential equations, (see Appendix, 
Arts. 102 and 103). 

61. The values to the fourth order are then obtained 
from those to the third by continuing the same process, and 

* 'En reflgchissant sur lea termes que doivent introdnire tontes les qnantit^s 

* pr6c^entes, on voit qu'il se pent gliaser des cosiniis de Tangle 6 dont nous avons 

* yn le dangereox effet d'amener dans la yaleor de u des arcs an lieu de leurs 
'cosinnsj de tels termes yiendront, par exemple, de la combinaison des cosinua 
'de (1 — m) avec des cosinns de m^ 

* Pour Iviter cet inconvenient qui 6terait & la solution pr6c6dento 

Tayantage de conyenir & un aussi grand nombre de r6yolutions qu*'on youdndt, 
*■ et la priyerait de la simplicity et de V uniyersalit6 si pr6cieuses en math^matiques, 
*il faut commencer par en chercher la cause. Or, on d^couyre facilement que 
' ces termes ne yiennent que de oe qu'on a suppose fixe Tapog^e du soleil, ce 

* qui n'est pas permis en toute rigueur, puisque quelque petite que soit sur cet astre 
' Taction de la lune, elle n'en est pas moins r^eUe et doit lui produire un mouyement 

* d'apog^ quoique trfes lent & la y^rit^'.— Clairaut, ThwrxA de la Lune, p. 55. 
2me Edition. 
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66 LUNAR THEORY. 

80 on to the fifth and higher orders; but the calculations are 
80 complex that the approximations have not been carried 
beyond the fifth order, and already the value of 6 in terms 
of t contains 128 periodical terms, without including those 
due to the disturbances produced by the planets. The coeffi- 
cients of these periodical terms are functions of 9n, e^ e'^ ~ , 

c, g^ ky and are themselves very complicated under their 
literal forms : that of the term whose argument is twice the 
difference of the longitude of the sun and moon, for instance, 
is itself composed of 46 terms, combinations of the preceding 
constants. 

See Pont^coulant, Syatime du Monckj torn, iv., p. 572. 
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CHAPTER V. 

NUMERICAL VALUES OP THE COEPPICIENTS. 

62. Having thus, from theory, obtained the form of the 
developments of the coordinates of the moon's position at 
any time, the next necessary step is the determination of the 
numerical values of the coefficients of the several terms. 

We here give three different methods which may be 
employed for that purpose, and these may, moreover, be 
combined according to circumstances. 

63. First Method. The values of the constants p, m^ a, 
)3, 7, \^ which enter into the arguments^ and of the additional 
ones a, A;, «, &c. which enter into the coefficients^ of the 
terms in the previous developments, may be obtained with 
great accuracy from observation. 

For this we must employ observations separated by very 
long intervals, such, for instance, as ancient and modern 
eclipses; also particular observations of the sun and moon 
made when the bodies occupy certain selected positions; 
and other observations of a special character. 

When the values of the elements have been so obtained, 
then the theoretical values of the coefficients may be com- 
puted by substitution in the analytical expressions. 

64. Second Method. Let the constants which enter into 
the arguments be determined as in the first method ; and let 
a large number of observations be made, from each of which 
« value of the true longitude, latitude, or parallax is obtained, 

I 
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58 LUNAR THEORY. 

together with the corresponding value of t reckoned from 
the fixed epoch when the mean longitude is zero. Let these 
corresponding values be substituted in the equations, each 
observation thus giving rise to a relation between the unknown 
constant coefficients, 

A very great number of equations having been thus 
obtained, let them, by the n^thod of least squares or some 
analogous process, be reduced to as many as there are 
coeflScients to be determined. The solution of these simple 
equations will give the required values. 

This method, however, would scarcely be practicable in 
a high order of approximation. For instance, in the fifth 
order, as stated in Art. (61), each observation would give 
rise to an equation containing 129 unknown quantities and 
the immense number of equations so obtained would have 
to be reduced to 129 equations of 130 terms each. 

65. Third Method. When the constants which enter into 
the arguments have been determined by the first method, we 
may obtain any one of the coefficients independently of all 
the others by the following process, provided the number of 
observations be very great. 

Let the form of the function be 

F=-4 + -Bsin^+ (7sin^ + &c, 
and let it be required to determine the constants -4, -B, (7, &c. 
separately ; ^, <^, &c., being known functions of the time. 

Let the results of a great number of observations corre- 
sponding to values 6^^ 6^, 0^^ &c., (f>^^ ^„ ^3, &c., beT,, F„ V^, 

&c. ; so that 

F^ = J+J?sin^j+(7sin0i + &c., 
V^ = A + Bs\ne^-\- (7sin</>2 + &c., 
V^ = A'^Bs\n0^+ Csin^3 + &c., 

F„ = ^ + i? sin ^,,+ Csin^^ + &c. 
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Now, n being very great^ we may assume that the sum of 
the positive values of each periodical term will be about 
counterbalanced by the sum of its negative values; and 
therefore, that if we add all the equations together these 
terms will disappear ; 

therefore A = — ^ ^ ^ — -_^^21L « 

n ' 

which determines the non-periodic part of the function. 

To determine B, Let the observations be divided into 
two sets separating the positive and negative values of sin^; 
then the other periodical terms, not having the same period, 
may be considered as cancelling themselves in adding up the 
terms of each set. Let there be r terms in the first set 

and s terms in the second, and let V\ F'', F' be the 

values of F corresponding to positive values of sin 9y which 
values we may assume to be uniformly distributed from 

fiin^ to siuTT, and therefore to be sin 8^, sin 28^, sinr.8^, 

where r . 85 = tt. 

And, again, let F^, F^^, F^^^, F,, be the values of F 

corresponding to the negative values of sin5 ; viz., sin (- a5), 
sin(-2A5), sin(— .9.A5), where 5.A5 = 7r. Then, 



F' =^ + 5sin85 + GAai/ +..., 
F'' = ul + J?sin2.S5+ (7sin^''+..., 

F' =^ + J?Binr.S0+ (7sin^*"+...; 
therefore 

F'+ F'+...+ F-=r.^+J?2,^(sin5) 

^tA^?^^: {And. hey, 

therefore 

T TtJo 

TT 



F^ =.4-J?sinA5 + (7sin<^^ +..., 
F,, = ^-J?sin2.A5+ (7sin^,,+..., 

V, =-4 — ji?sins.A5+ (7sin^, +...; 

therefore 

F, + F,,+...+ F=5.^-J?V(8me) 

= 5^--^V(sIn5.A5); 
therefore 



— ' " =A — I sin p . 



de 
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therefore B-- 



1\ r 8 J^ 

and in a similar manner may each of the coefficients be 
independently determined.* 

66. When the periods of two of the terms differ but 
slightly — for instance if and ^ go through their periodic 
variations very nearly in the same time,^— the method could 
not then with safety be applied ; for, since the same values 
of and <f> would very nearly recur together during a longer 
time than that through which the observations would ex- 
tend, the two terms would be so blended in the value of 
Fthat they would enter nearly as one term — the difference 
between and ^ would be very nearly the same at the 
end as at the begmning of the series of observations. 

* If r and s acp not Bofficlently great to allow ns to subetitnte f^'^^'anBdO for 
£«'"' sin • ^^f we must proceed as follows ; 

F' + F" + ,+ V''=rA + B (sin^O + sm2de + + sinrae) 

V'+V"+. + V' . B l~it9' , 

Similaily, -^ ii-j i = A--^^l^if^j; 

therefore 

nr f V'+V»»+ ^yr V^+V„+ + FA 

"*{^-if^*ir ' - ■ > 
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67. Let US suppose the periods to be actuallj Identical^ 
so that ^ s d + a, a being some constant angle ; then 

BemO + Cem^ 

may be written (J5+ (7cosa) sin^+ Osina cos^, 

or F=-4 + (J5+(7cosa)sm<?+08inacos(?+ * 

If now we divide the observations, as before, into two 
sets, corresponding to the positive and negative values of 
sind, the terms involving cosd will disappear in the sum- 
mation of each set ; and, following the process of the method, 
we shall find the value of 

jB + C7 cos a = Jf suppose. 

Dividing again into two sets corresponding to the positive 
and negative values of cos 9, the terms in sin d will be can- 
celled, and the same process will give the value of 

(7 sin a =s -N" suppose. 

Treating the observations in the same way with respect 
to the angle ^, we get two results, 

(7+£cosa = 2f', 

-^sina-JV^'; 

from these four equations we easily get 

Mj Nj M'j N' are connected by the equation of condition, 

68. When the periods of and ^ are nearly, but not 
exactly, the same, this equation of condition will not hold, 
and the preceding ^values of B and C would not be exactly 
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correct ; but yet they would be very approximatCi especiafly 
if the mean between the two values of B be taken. 

Or, we may, after having taken one of these slightly 
erroneous values for J?, make a further correction by estab- 
lishing as it were a counterbalancing error in the value of 
C, Let R be the value so found for J?; then, from the 
V of each of the observations subtract the value B^ sin 6^ 
the result Z7will be very nearly equal to -4 + (7 sin + &c., 
and from the n equations 

l7^ = -4-f>Osin0,+ 

J7, = ^+(7sin^,+ 

?7^=^+asin^„+ 

a value (7 of G will be obtained, by the rule, which will be 
very approximate, and, at the same time, agree better with 
jB' in satisfying the equations than G itself would do. 

69. When two terms whose periods are nearly equal do 
occur, it is plain, by examining the values of M and M\ that 
the errors which would be committed by following the rule, 
without taking account of this peculiarity, would be the taking 
B+ G cosa aniG+B cos a for B and C respectively. 
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CHAPTER VL 



PHYSICAL INTBBPBETATION. 



70. The solution of the problem which is the object of 
the Lunar Theory may now be considered as effected ; that 
is, we have obtained equations which enable us to assign 
the moon's position in the heavens at any given time to 
the second order of approximation ; we have explained how 
the numerical values of the coeflBcients in these equations 
may be determined from observation; and we have, more- 
over, shewn how to proceed in order to obtain a higher 
approximation.* 

It will, however, be interesting to discuss the results we 
have arrived at : — to see whether they will enable us to form 
some idea of the nathre of the moon's complex motion ; also 
how far they will explain those inequalities or departures 
from uniform circular motion which ancient astronomers had 
observed ; but which, until the time of Newton, were so 
many unconnected phenomena; or, at least, had only such 
arbitrary connexions as the astronomers chose to assign, by 
grafting one eccentric or epicycle on another as each newly 
discovered inequality seemed to render it necessary. 

It is true that our expressions, composed of periodic 
terms, are nothing more than translations into analytical 

♦ The means of taking into account the ellipsoidal figure of the earth and the 
disturbances produced by the planets, are too complex to form part of an intro- 
ductory treatise. For information on these points reference may be made to 
Airy's Figure of the Earth, Pont^coulant's Systeme du Monde, vol. iv. 
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64 LUNAR THEORY. 

language of the epicycles of the ancient;* but they are 
evolved directly from the fundamental laws of force and 
motion, and as many new terms as we please may be ob- 
tained by carrying on the same process ; whereas the epicycles 
of Hipparchus and his followers were the result of numerous 
and laborious observations and comparisons of observations ; 
each epicycle being introduced to correct its predecessor 
when this one was found inadequate to give the position 
of the body at all times : just as with us, the terms of the 
second order correct the rough results given by those of 
the first;* the terms of the third order correct those of the 
second, and so on. But it is impossible to conceive that 
observation alone could have detected all those minute 
irregularities which theory makes known to us in the terms 
of the third and higher orders, eyen supposing our instruments 
far more perfect than they are ; and it will always be a 
subject of admiration and surprise, that Tycho, Kepler, and 
their predecessors should have been able to feel their way 
so far among the Lunar inecjualities, with the means of 
observation they possessed. 

LONGITUDE OP THE MOON. 

71. We shall first discuss the expression for the moon^a 
longitude, as found Art. (55), 

=:pt + 2e sin {cpt - a) + {e' sm2 {cpt- a) 
+ ^me sin {(2 - 2wi - c)p* - 2)8 + a} 
+ V^* »"» {(2 - ^m)pt - 2)8} 

- ^me' Bin{mpt-^ fi - ^) 

— lA*«wl2(5p^-7)* 
The mean value of ^ is p^; and in order to^ judge of the 
effect of any of the small terms, we may consider them 

♦ See Whewell'a History of the Inductive Sciences, 
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ELLIPTIC INEQUALITY. 65 

one at a time as a correction on this mean valae pty or 
we may select a combination of two or more to form this 
correction. 

We shall have instances of combination in the elliptic 
inequality and the evection^ Arts. (73) and (77) ; but in the 
remaining inequalities each term of the expression will form 
a correction to be considered by it&ielf. 

72. Neglecting all the periodical terms, we have 

de 

which indicates uniform angular velocity; and as, to the 
same order, the value of u is constant, the two together 
indicate that the moon moves uniformly in a circle, the 

period of a revolution being — , which is, therefore, the 

expression for a mean sidereal month, or about 27 J days.* 
The value of^ is, according to Art. (54), given by 

and 2A m is due to the disturbing action of the sun, we see 
that the mean angular velocity is less, and therefore the mean 
periodic time greater than if there were no disturbance. 

Elliptic Inequality or Equation of the Centre. 

73. We shall next consider the effect of the first three 
terms together: the effect of the second alone, as a correction 
oipt^ will be discussed in the Historical Chapter, Art. (115). 



♦ The accurate value was 27d. 7h. 43m, ll'261s. in the year 1801. See 
Art. (108). 

K 
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66 LUNAR THEORT. 

=zpt + 2e sin {cpt - a) + Je* sin 2 (cpt - a), 

which may be written 

(?=p«+2esm[j?«--{a+(l-c)p«}]+fe*Bin2[p«-{a+(l-c)j)<}]. 

But the connexion between the longtitnde and the time in 
an ellipse described about a centre of force in the focus is. 
Art (13), to the second order of small quantities : 

^ = n« + 2^8in(ne-a') + f(g*8in2(n^-a'), 

where n is the mean motion, e the eccentricity, and a' the 
longitude of the apse.* 

Hence, the terms we are now considering indicate motion 
in an ellipse ; the mean motion being p, the eccentricity «, 
and the longitude of the apse a + (l — c)^^; that b, the 
apse is not stationary but has a progressive motion in longi- 
tude, uniform, and equal to (1 — c)p. 

74. The two terms 2e sin (cpt — a) + f e" sin 2 (cpt — a) con- 
stitute the elliptic inequality^ and their effect may be further 
illustrated by means of a diagram. 

Let the full line AMB re- 
present the moon's orbit about 
the earth J5', when the time t 
commences, that is, when the 
moon's mean place is in the prime 
radius -ffT, from which the lon- 
gitudes are reckoned. 

The angle TZ4, the longi- 
tude of the apse, is then a. At the time e, when the moon's 
mean longitude > r^JIf=^^, the apse line will have moved 




♦ The epoch c which appears in the expression of Art. (18) is here omitted; 
a proper aflsumption for the origin of t, as explained in Art. (86), enabling q^ 
to avoid the e. 
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MOTION OP THE APSE. 67 

in the same direction through the angle AEA =£ (1 — c) TEM^ 
and the orbit will have taken the position indicated by the 
dotted ellipse. The true place of the moon in this orbit, 
so far as these two terms are concerned, will be m, where 

MEm = 26 sin [cpt — a) + fe* sin 2 [iq^t — a) 

=3e sin^'jEM + Je" sin 2AEM . 

= 26 mi A EM {\ + fe Qo^AEM) ; 

which, since e is about ^, is positive from perigee to apogee, 
and therefore the true place before the mean; and the 
contrary from apogee to perigee: at the apses the places 
will coincide. 

76. The angular velocity of the apse is (1 — c)p, or, if 
for c we put the value found in Art. (52), the velocity will 
be fm'jp. Hence, while the moon describes 360", the apse 
should describe |m'.360* = If* nearly, m being about -j^. 

But Hipparchus had found, and all modem observations 
confirm his result, that the motion of the apse is about S"* in 
each revoludon of the moon. See Art. (118). 

This difference arises from our value of c not being repre- 
sented with sufficient accuracy by 1 — |m*. 

Newton himself was aware of this apparent discrepancy 
between his theory and observation ; and we are led, by his 
own expressions (Scholium to Prop. 35, lib. iii. in the first 
edition of the iWnctjpta), to conclude that he had got over 
the difficulty* This is rendered highly probable when we 
consider that he had solved a somewhat similar problem 
in the case of the node ; but he has nowhere given a state- 
ment of his method : and Clairaut, to whom we are indebted 
for the solution, was on the point of publishing a new hy- 
pothesis of the laws of attraction, in order to account for 
it, when it occurred to him to carry the approximations to 
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68 LUNAR THEORY. 

the third order, and he there found the next term in the 
value of c nearly as considerable as the one already obtained. 
See Appendix. The value is 

/. (1 - c) 360° = (14 /^^) (value found previously) 
= 2|** nearly, 

thus reconciling theory and observation, and removing what 
had proved a great stumbling-block in the way of all as- 
tronomers.* 

When the value of c is carried to higher orders of ap- 
proximation, the most perfect agreement is obtained. 

The motion of the apse line is considered by Newton 
in his Brinciptaj lib. I., Prop. 66, Cor. 7. 

Evection. 

76. The next term + J^wi6sin{(2- 2m-c)p«-2^ + a} 
in the value of 6 has been named the Evection. We shall 
consider its effect in two different ways. 

First, by itself, as forming a correction on pt, 

=zpt + ^me sm {(2 - 2m - c)^« - 2^ + a\. 
Let D = pt s= moon's mean longitude at time t^ 

0= mpt+P c= sun's ...•••.. ..•.•»».*,.••., 

a' == (1 - c) pt + a smean longitude of apse .*......., 

then 



6 =pt + ^me sin [2 {pt - {mpt + fi)} - {pt -{l-c)pt + a}] 
=i>f + ^me sin (2 (J - o) - (J - a')}. 

* See Dr. Whewell's Bridgwater Treati$6, 
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EVECTION. 69 

The effect of this term will therefore be as follows : 
In syzygies 

d^pt-^ ^me sin (D - a') ; 

or the true place of the moon will be before or behind the 
mean, according as the moon, at the same time, is between 
apogee and perigee or between perigee and apogee. 
In quadratures 

=zpt + ^me sin (J — a'), 

and the circumstances will be exactly reversed. 

In both cases, the correction will vanish when the apse 
happens to be in syiaygy or quadrature with regard to the 
sun at the same time as the moon. 

In intermediate positions, the nature of the correction 
is more complex, but it will always vanish when the sun 
is at the middle point between the moon and the apse, 
or when distant SO"" or 180'' from thfit point; for if 

J + a' 
G = ^ r.90 , where r = 0, ± 1, or 2, 

sin[2(>-O)-(l)-a')] = sin(l) + a'-20) 
=rsinr.l80'' 
=0. 

77. The other and more usual method of considering 
the effect of this term is in combination with the two terms 
of the elliptic inequality, as follows : 

To determine the change in the position of the apse and in 
the eccentricity of the moovUs orbit produced by the evection. 

Taking the elliptic inequality and the evection together^ 
we have 

=ipt + 2e sin (cp« — a) + f €* sin2 {ppt - a) 

+ ^me sin {(2 - 2m- c);?«- 2i8 + a}. 
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70 LUNAR THEORY. 

* 

Let ci be the longitude of the apse at time t on sappositioii 
of uniform progression, 

• sun ;* 

whence a' = (1 — c)pt + a, 

. 0= mpt +)8. 
And the above may be written 

sxpt + 2e sin [cpt — a) + Je* sin2 {cpt — a) 

+ ^me sin {g^- a + 2 (a' - ©)} ; 
and the second and fourth terms may be combined into one, 

2-Bsin(cpe-a + S), 
if we assume jEcosS = e + ^me cos2 (a' — o)i 
and J&sinS= y^7nesin2(a'- o); 

whence tan 8 and ^may be found; and approximatelji 
E= e {1 + ^m cos2 (a' - o)}, 
S=^msin2(a'-0). 

The term fe' sin2(cp<— a) may also, therefore, to the 
second order, be expressed by 

fj&'sin2(<ye-a + S}, 

and the longitude becomes) 

6=^pt + 2jFsin ((ye - a + 8) + |-B* sin2 {cpt-a+ 8), 

or 0^pt + 2JE?sin {pt - a' + 8) + f JE?" sm2 (p< - a' + 8) j 

but the last two terms constitute elliptic inequality in an orbit 
whose eccentricity ia E and longitude of the apse a' — 8 ; 
therefore the evection, taken in conjunction with elliptic 
inequality, has the effect of rendering the eccentricity of the 
moon*8 orbit variable, increasing it by ^me when the apse- 
line is in syzygy, and diminishing it by the same quantity 
when the apse-line is in quadrature ; the general expression 
for the increment being ^me cos 2 (a' — o). 
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And another effect of this term is, to diminish the longitude 
of the apse, calculated on the supposition of its uniform 
progres^on, by the quantity 8 = y^wisin2 (a'- O); so that 
the apse is behind its mean place when in the first and third 
quadrant in advance of the sun, and before its mean place 
in the second and fourth. 

The cycle of these changes will evidently be completed 
in the period of half a revolution of the sun with respect 
to th^ apse, or in about -j^ of a year.* 

78. The period of the eveotion itself, considered inde- 
pendently of its effect on the orbit, is the time in which 
the argument (2 — 2w — c)pt — 2)8 + a will increase by 27r. 

Therefore period of evection 

^ 27r , _ mean sidereal month 

"" (2-2wi-c)jp "" 2-2m-c 

mean sidereal month 



=a^,.«riy. 



l-2m + |m' 
== 31^ days, nearly.f 

Newton has considered the evection, so far as it arises 
from the central disturbing force, in Prop. 66, Cor. 9, of 
the Principia. 

Variatian. 
79. To explain the physical meaning of the term 
V^' sin {{^- 2m) pe - 2i8}, 

* The change of eccentricity and the yariatipn in the motion of the apse 
follow the same law as the abscissa and ordinate of an ellipse referred to 
its centre: for if £ - e = a; and d = y, then 

t The accurate Talue is 31 * 8119 days. 
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72 LUNAR THEOEY. 

in the expression Jbr the moorCs longitude^ 

e ^pt + V^* sin {(2 - 2m) pt - 2)8}. 

Let }) represent the moon's mean longitude at time t^ 

O • san*s •••• 

therefore J —pt^ 

Q^mpt + fi; 

and the valae of becomes 

^=p<+V^*8in2(l)-0), 

which shews that from syzjgy to quadrature the moon's 
true place is before the mean, and behind it from quad- 
rature to sjzygy; the maximum difference being ^m* in 
the octants. 

The angular velocity of the moon, so far as this term 
is concerned, is 

■^ =|> + V (1 - ^) ^'i' C082 (J - O), 

■ep {1 + ^m' cos2 (J - o)}, nearly, 

which exceeds the mean angular Telocity j? at syzygies, is 
equal to it in the octants, and less in the quadratures. 

This inequality has been called the Variation^ its period 
is the time in which the argument (2 — 2m) pe — 2)3 will 
increase by 27r; 

. , ^ . ^. 27r mean synodlcal month 
.-. penod of vanation = -(2^2^ —2 

= 14| days, nearly.* 

80. The quantity ^m* Is only the first term of an endless 
series which constitutes the coefficient of the variation, the 
other terms being obtained by carrying the approximation 

* The accurate value is 14*765294 days. 
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to a higher order. It is then found that the next term in 
the coefficient is f |w', which is about ^^ of the first term ; 
and as there are several other important terms, it is only 
by carrying the approximation to a higher order (the 5th 
at least) that the value of this coefficient can be obtained 
with sufficient accuracy from theory. In fact, ym" would 
give a coefficient of 26' 27'' only ; whereas the accurate value 
is found to be 39' 30". 

The same remark applies also to the coefficients of all 
the other terms. 

81. As far as terms of the second order, the coefficient 
of the variation is independent of e the eccentricity, and k 
the inclination of the orbit. It would therefore be the same 
in an orbit originally circular, whose plane coincided with 
the plane of the ecliptic: it is thus that Newton has con- 
sidered it. Princip. prop. 66, cor. 3, 4, and 5. 

Anmial Equation. 

82. To explain the physical meaning of the term 

— 3me' sin [mpt + ^8 - f ) 
in the expression for the moon's longitude, 
^ =pt — 3iwe' sin {mpt + iS - f ), 

=pt — 8wte' sin (longitude of sun - longitude of sun's perigee), 

^pt — 3me* sin (sun's anomaly). 

Hence, while the sun moves from his perigee to hi& 
apogee, the true place of the moon will be behind the mean ; 
and from apogee to perigee, before it. The period being an 
anomalistic year, the effect is called Annual Equation. 

Differentiating we get 

-J —p [\ — 3mV cos (sun's anomaly)}. 
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74 LUNAR THEORY. 

Hence, so far as this inequality is concerned, the moon^s 
angular velocity is least when the sun is in perigee, that 
is at present about the Ist of Januaryjand greatest when 
the sun is in apogee, or about the 1st ctfjuly. 

The annual equation is, to this order, independent of the ^ 
eccentricity and inclination of the moon's orbit, and therefore, 
like the variation, would be the same in an orbit originally 
circular. Vide Newton, Principia^ prop. 66, cor. 6. 

Eeduction, 

83. Before considering the effect of the term 

A* 
--sin2(gp«-7), 

which, as we shall see Art. (89), is very nearly equal to 
the difference between the longitude in the orbit and the 
longitude in the ecliptic, it will be convenient to examine 
the expression ^for the latitude of the moon, and to see how 
the motion of the node is connected with the value of g. 

LATITUDE OF THE MOON. 

84. The expression found for the tangent of the latitude,* 
Art. (51), is 

s = A;sin(^^-7)4|w2Asin{(2-2m-^)5-2^ + 7}. 

If we reject all small terms, we have 

5 = 0, 

or the orbit of the moon coinciding with the ecliptic, which 
is a first rough approximation to its true position. 



* This expression for the tangent of the latitude is more convenient than 
that which gives it in terms of the mean longitude, Art. (57), on account of 
the less number of terms involved. See Pontecoulant, vol. iv. p. G30. 
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LATITUDE OF THE MOON. 75 

85. Taking the first term of the expansion 
. s = Asin(5'^ — 7), 
we may write it 

Let TNm be the ecliptic, N the 
moon's node when her true longitude 
is zero, and let M be the position 
of the moon at time ty m her place 
referred to the ecliptic ; 

therefore T-A7'«=7, Twi = ^, tB,nMm^s, 

Take Nlf' = {g-1)0 in a retrograde direction^ and join 
MN^ by an arc of great circle ; 

then tan Mm = tan MWm sin Wm , 

or 5 = tan MN'm sin [^ - {y - (^ - 1) 0]] ; 

which, compared with the value of s given above, shews that 
Mif^m = tan**A is constant, and therefore the term k sin {g0—y) 
indicates that the moon moves in an orbit inclined at an 
angle tan'^i to the ecliptic, and whose node regredes along 

the ecliptic with the velocity (S' - 1) -ii > ^^ ^^^^ ^ mean 

velocity (5^— 1)^?. 

86. Hence the period of a revolution of the nodes 
_ 27r _^ one sidereal month 

but, from Art. (51), the value of ^r = l + Jm* ; 

^, ^ • 1 /• 1 .- n 3 one sidereal month 
therefore period of revolution of nodes = s— = 

= 6511 days, nearly. 

This will, for the same reason as in the case of the apse, 

Art. (75), be modified when we carry the approximation 
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to a higher degree ; this value of g Is, hawever, much more 
accurate than the corresponding value of c, for the third 
term of g is small ; the value to the third order being (see 
Appendix, Art. 102) 

, ,, . 1 /. 1 .. /. .1 1 one sidereal' month 
and the period of revolution of the nodes = — » .,, — 

= 6705 days, nearly. 

This is not far from the accurate value as given by 
observation, and when the approximation to the value of g 
is carried to a higher order, the agreement is nearly perfect. 

The true value is 6793-39 daySj that Is about ISyrs. 7 mo. 

Evection in Latitude. 

87. To explain the variation of the inclination and the 
irregularity in the motion of the node expressed by the term 

+ ^mk sin {(2 - 2m - ^) ^ - 2^8 + 7}. 

This term, as a correction on the preceding, is analogoug 
to the evection. as a correction on the elliptic inequality. 
Taking the two terms together, 

8 = k sin (gd '-y)+ ^mk sin {(2 -2171- g)0 - 2fi -{-ylf 

Let J = longitude of moon = 0y 

0= sun =^md + fii 

Q= node =7 — (5^-!)^; 

therefore « = A; sin (}) - g^ ) + ^mk sin {J - g^ - 2 (O - Q,)}. 

Now these two terms may be combined into one, 

« = Z'sin(l)- Q -A), 

if 5:'cosA = 4 + fmAcos2(0- Q), 

jS: sin A = Imk 8in2 (O - q), 
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whence A and K may be found ; and, approximately, 

ir=i{l + f7/*cos2(o-a)}, 
A = |?nsm2(0- Q), 
but the equation 

« = jr8in(})- Q- A) 

represents motion in an orbit inclined at an ai^le tan'^JT to* 
the ecliptic, and the longitude of whose node is Q + A. 

This term has therefore the following effects : 

Ist. The inclination of the moon's orbit is variable, its 
tangent increases by \mk when the nodes are in syzygies, 
relatively to the sun, and decreases by the same quantity 
when they are in quadrature; the general expression for 
the increase being \mk C082(o — Q^). 

2nd. The longitude of the node, calculated on^ supposition 
of a uniform regression, is increased by A = |m sin2 (o — ^.), 
so that the node is before its mean place while in^ the first 
or third quadrant in front of the sun, and behind it in the 
second and fourth. Prind'pia^ book ill., props. 33 and. 35; 

The cycle of these changes will be completed in the 
period of half a revolution of the sun with respect to lie 
node, that is, in 173*21 days^ not quite hsdf-a-year. 

88i The tangent of the latitude has here been obtained ^ 
if we wish to* have the latitude itself it will be given by the 
formula 

latitude = s - J^'* + \i- &c. 
which, to the degree of approximation adopted, will clearly 
be the same as «. 

lieducticm. 

89. We may now consider the term which we had 
neglected (Art. 83) in the expression for the longitude, namely, 
-iA*8in2(52>^-7). 
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78 LUNAR THEORt. 

Let N be the position of the node when th« moon^s 
longitude is 0j M the place of the 
moon, m the place referred to the 
ecliptic. 

Therefore Ttw = 6^ 

Nm=g0 — y^ 
tm N=k. 

The right-angled spherical triangle NMm gives 
tan -Win 




therefore 



cos^=^ tt^^, 

cos N tan NM - tan Nm 



1 + cos^ tan NM^ tan Nm ' 



^^ 2 " Bm(NM+Nm)' 

or, since both N and NM— Nm are small, 

tan'iV- NM-- Nm . , , 

therefore NM-Nm^^k^ sm2(g0—y)=lk^ Bm2(gpt'-y\ nearly. 

Hence this term, which is called the reduction^ is approxi- 
mately the difference between the longitude in the orbit and 
the longitude in the ecliptic. 

RADIUS VECTOR. 

90. To eocplain the physical meaning of the terms in the 
valtie of u,. 

We shall, for the explanation, make use of the formula 
which gives the value of u in terms of the true longitude, 
Art. (52> 
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TERMS IN RADIUS VECTOR. 79 

NeglectiDg the periodical terms, we have for the mean 
value 

The term — ^wi*, which is a consequence of the disturbing 
effect of the sun, shews that the mean value of the moon's 
radius vector, and therefore the orbit itself, is larger than if 
there were no disturbance. 

Elliptic Inequality. 

91. To explain the effect of the term of the first order ^ 

w = a {1 + e cos (cfl — a)}, 

= a[l + ecose-{a + Cl-cJ^}]. 

This is the elliptic inequality^ and indicates motion in an 
ellipse whose eccentricity is e and longitude of the apse 
a + (l — c)fl; and the same conclusion Is drawn with respect 
to the motion of the apse as In Art. (73). 

Evection, 

92. To explain the physical meaning of the term 

V'*^ cos {(2 - 2wi - c) - 2/3 + a}. 

This, as in the case of the corresponding term in the longi- 
tude, is best considered in connection with the elliptic 
inequality, and exactly the same results will follow. 

Thus, calling J, O, and a' the true longitudes of the moon, 
sun, and apse (the last calculated on supposition of uniform 
motion), these two terms may be written, 

i4 = a[l + ecos(])-a')- V^^<^^s{})--a' + 2(a'-0)}] 
= a[l + jE;cos(])-a' + S)]; 
where jE^ cos S = e + ^ me cos 2 (a' — 0), 

jE? sin S = ^me sln2 (a - 0). 
These are identical with the equations of Art (77). 
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80 LUNAR THEORY. 

Variation. 

93. To explain the effect of the term m^a cos {(2-2m) B- 2/8} , 

tt = a[l + w*co8{(2-. 2m)0-2^]] 
« a [1 + w* C082 (J - o)]. 

As far as this term is concerned, the moon's orbit would 
be an oval having its longest diameter in quadratures and 
least in syzygies. Principia^ lib. I. prop. 66, cor. 4. 

The ratio of the axes of the oval orbit will be 

^j 1 = IS nearly, m being "0748. 

See iVtnajpia, lib. iii. prop. 28. 

Reduction. 

94. The last important periodical term in the value of u is 

- — cos2(5re-7). 

This term expresses approximately the variation in the 
difference between the values of u in the orbit and in the 
ecliptic. 

For if u^ be the reciprocal of the value of the radius 
vector in the orbit, 

tt| = te cos (latitude), 

therefore w — w^ = Jt^* «= \ak^ sin* {gd — 7) 

= iaA* - lak^ cos 2 {gO - 7) 

«= const. — \ak* cos 2 (^gO — 7). 

Taking this result in connexion with that of Art. (89) 
we see that the reduction^ so far as periodical terms to the 
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second order are concerned, is simply a geometrical con- 
sequencejof the inclination of the orbit ; and that, if we 
measured the longitude and the radius vector along the 
orbit instead of taking their projection on the plane of the 
ecliptic, these periodical terms would not appear. 

95* There are no other terms of the second order in the 
value of u. The annual equation^ which, in the longitude, 
is of the second order, is only of the third order in the 
radius vector. 

Periodic time of the Moon. 

96. We have seen. Art. (72), that the periodic time of 
the moon is greater than if there were no disturbing force ; 
but this refers to the mean periodic time established on an 
interval of a great number of years, so that the circular 
functions in the expression are then extremely small com- 
pared with the quantity pt which has uniformly increased. 

When, however, we consider only a few revolutions, these 
terms may not all be neglected. The elliptic inequality and 
the evection go through their values in about a month, the 
variation and redu>ction in about half-^a-*month ; their effects, 
therefore, on the length of the period can scarcely be con-^ 
sidered, as they will increase one portion and then decrease 
another of the same month. 

But the annttal equation takes one year to go through its 
cycle, and, during this time, the moon has described thirteen 
revolutions; hence, fluctuations may, and, as we shall now 
shew, do take place in the lengths of the sidereal months 
during the year. 

We have, considering only the annual equation, Art. (82), 

pt=^e-^ Sme sin (tw^ + /3 - ?)• 

M 
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82 LUNAR THEORY. 

Let T be the length of the period: then, when is 
increased by 27r, t becomes ^ + jT; 

therefore ^ (< + T) = 27r + 5 + %me! sin (2m7r + m^ + /S - ?), 
whence ^ r= 27r + 6?we' sin twtt cos (ww + m^ + i8 — 5 

therefore T— mean period + sinrnzr cos(o - f), 

where O =7n^ + i8 + 7n7r = sun's longitude at the beginning 
of the month + twtt 
= sun's longitude at the middle of the month. 
Hence Twill be longest when G - ?=0, 
and shortest when O — f = tt ; 

or T will be longest when the sun at the middle of the month 
is in perigee, and shortest when in apogee ; but, at present, 
the sun is in perigee about the 1st of January, and in apogee 
about the 1st of July; therefore, owing to annual equation^ 
the winter months will be longer than the summer months, 
the difference between a sidereal month in January and July, 
from this cause, being about 20' minutes. 

97. All the inequalities or equations, which our expres- 
sions contain, have thus received a physical interpretation. 
They were the only ones known before Newton had estab- 
lished his theory, but the necessity for such corrections was 
fully recognized, and the values of the coefficients had already 
been pretty accurately determined ; still, with the exception 
of the reduction, which is geometrically necessary, they were 
corrections empirically made, and it was scarcely to be ex- 
pected that any but the larger inequalities, viz. those of the 
first and second orders which we have here discussed, could 
• be detected by observation: we find, however, that three 
others, have, since Newton's time, been indicated by obser- 
vation before theory had explained their cause. These are — 
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the secular acceleration^ discovered by Halley ; an inequality, 
found by Mayer, in the longitude of the moon, and of which 
the longitude of the ascending node is the argument ; and 
finally an inequality discovered by Btirg, which has only of 
late years obtained a solution. For a further account of 
these, as also of some other inequalities which theory has 
made known, see Appendix, Arts. (108) to (112). 
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CHAPTER VII. 



APPENDIX. 



In this chapter will be found collected a few propositions 
intimately connected with the results or the processes of the 
Lunar Theory as explained in the previous pages. Befer- 
ence has been made to some of them in the course of the 
work, and the interest and importance of the others are 
sufficient to justify their introduction here. 

98. The moon is retained in her orbit hy the force of 
gravity^ that is; by the same force which acts on bodies at 
the surface of the earth. 

The proof of this is merely a numerical verification ; the 
data required from observation are, 

the space fallen through from rest in 1" by bodies at the 

earth's surface = 161 feet, 

the radius of the earth = 4000 miles, 

the periodic time of the moon = 27j days, 

the dbtaace of the moon from the earth's centre=60x4000 miles. 

The force of the earth's attraction oc -rr. — r^ . Therefore, the 

(dist.)' ' 

space fallen through in 1" at distance of moon by a body 

1 fi*i 
moving from rest under the earth's action = -^ feet 

= -00447 feet. 
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But the moon in one second describes an angle . ^ =ca, 
during which the approach to the earth 

= 60 X 4000 X 5280 (vers, w) feet 
60 X 4000 X 5280 x 27r' 



(27J/.(24)^(60/ 
= -00448 feet. 



feet 



( Therefore, the space through which the moon is deflected 
in one second from her straight path, is just the quantity 
through which she would fall towards the earth, supposing 
her to be subject to the earth's attraction; and we may, 
therefore, conclude that she is retained in her orbit by the 
force of gravity. 

When first Newton, in 1666, attempted to verify this 
result, he found a difference between the two values equal 
to one-sixth of the less: the reason of his failure was the 
incorrect measures of the earth, which he made use of in 
his computation; and it was not till about 16 years later 
that he was led to the true result, by using the more 
correct value of the earth's radius obtained by Picart. 

Principiaj lib. III., prop. 4. 

99. The moon's orbit is everywhere concave to the sun. 

Let 8y Ey and M be the centres of the sun, earth, and 
moon. We ought first to apply 
to each body forces equal and 
opposite to those which act on 
the sun in order to bring him to rest. These forces are, 
however, so small that we may neglect them, and we shall 
consider the moon to be moving about the sun fixed, and to 
be disturbed by the earth s^lone. 
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The forces on M are, therefore, -oijts * 






.in Jf5, 
.in ME. 



This last must be resolved Into two, one in MS^ the other 
perpendicular to it. 

Therefore, the whole central force on the moon in M8 



m 



E 



and the proposition will be proved if we shew that this force 
is always positive. Now, Art. (25), 



periodic time of sun = 

and 

therefore 



2ir.8G^ 27r.SM^ 



m^ 



.moon 



27r.EM^ 



nearly. 



therefore 



therefore 



"(E+My 

m' _ E+M . 
^jy^s-T*^-;^]^ nearly; 

E 



m 
SM 

SM'' 



8 > tJ-? ^Jj^» i 

. SM ^ E 



EM' EM* 



EM' 



E 



SM^'EM*''^'"''^'''^'' 

but this is the value of the central force corresponding to 
cosilf=— 1, and is therefore its least value. Hence the 
central force always tends to the sun, and the path is always 
concave. 

At new moon the force with which the moon tends to the 
sun is, therefore, greater than that with which she tends to 
the earth: the earth being itself in motion in the same 
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CENTRAL DISTURBING FORCE. 87 

direction, and, at that instant, with greater velocity, will 
easily explain how, notwithstanding this, the moon still 
revolves about it.* 

Central and Tangential Disturbing Forces, 

100. We have hitherto considered the effects of the 

central and tangential disturbing forces in combination ; but 

it will be interesting to determine to which of them the 
several inequalities principally owe their existence. 

(I) To determine the effect of the central disturbing force. ' 

Make r=0; 

therefore ^^^^■^-F^? = ^' 

p 

or substituting for j^^ from Art. (52), 



f 1 - W - i^' + I ^^^ ^03 ^^1 - f ^' cos (2 - 27n) d^ 
^'!^ + M=a -f|m'ecos(2-2w-c)^j 

"^^ I +iA*cos2^^,; 

therefore, 

a - p'- Jm'4e cos(ce- a) +im''cos[(2 - 2m)^-2/3} 
u=a\ IH- -^^me cos {(2 - 2?n - c) ^ - 2^8 + a} 

I -p' cos 2 (^^-7). 

If we compare this with the value of u found Art. (52), we 
see that the elliptic inequality and the reduction are due 
to the central or radial force, as also one-half of the varia- 
tion and three-tenths of the evection. ^ ^ 

It would perhaps be proper to separate the absolute cen- 
tral force from the central disturbing force; the terms due 
to the latter are those which contain m ; therefore, the elliptic 

* See the Author's Astronomy^ p. 233. 
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inequality and the reduction are the effects of the former, 
except that in the elliptic inequality the introduction of c, or 
the motion of the apse, is due to the disturbing force. 

(2) To determine ike effect of the tangential disturbing force. 

Let the central disturbing force be zero ; 

P II 
then jj-^ = p (1 — |s*) = a, neglecting the inclination, and, 

substituting for ^, . £ and for 2 (^ + «) f^. d0 their 

values from Art. (52), the equation in u becomes 

d*u 

^^+w = a{l-f7w*cos(2-2w)^j + ^w*ecos(2-2»i-c)^J; 

, _ fH-ecos(cfl-a) + im'cos{(2-2in)^-2^} 

wnencew-aj ^.|^^geos{(2-2m-c) ^-2/3 + a}. 

We have here the remaining half of the variation and the 
remaining seven-tenths of the evection. as the effects of the 
tangential disturbance. Also c = l, or, to the second order, 
the tangential force has no effect on the motion of the apse. 

101. The separate effects of the central and tangential 
forces in producing the inequalities in the longitude may be 
traced in a similar manner. We shall leave this as an 
exercise for the student. The annual equation will be found 
to be due to the central force. The other inequalities will 
be divided as in the radius vector; except as regards the 
variation, of which four-elevenths are due to the central 
and the remaining seven-elevenths to the tangential force. 

To calculate the value of g to the third order. 

102. We must here make use of the results which the 
approximations to the second order have furnished; but as 
the value of g is determined by that term of the differential 
equation whose argument is gO^-yy we need only consider 
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those terms which by their combinations will lead to It with- 
out rising to a higher order than the fourth. 
From Arts. (51), (31), we obtain 

8=^k siugd^ + imk sin (2 - 2m -g) 0^^ 

^^V~="-|»«^{l+cos(2-27n)^J, 

T 
^8 = -Xsin(2-2m)^,. 

From these we obtain 

Ps- 8 

-^-^ =.••- f »i*^ (1 - y^m) sin^r^^ 

T ds 
'Wi^^db =—"f ^"*sin (2 -2w) d^ {cos^r^^+fwi cos(2-2w-5r) ^J 

=...--^^WAsin5r^j, 

d*8 

~- + 5 = 0, to the second order : 

therefore 2 f ^ + « J It,— » ^^ = 0, to the fourth order. 

Therefore the equation in 5 becomes, so far as these terms 
are concerned, 

assume « = A sin^^^j, 

therefore A (1 - /) = - f wi'A + ^^ m'*, 

<7 = l + |m'-^5m'. 

To find the value of c to the third order. 

103. Proceeding in a similar manner from the equations. 
Arts. (52), (31), we get 

tt=a{l+ecosc^i+w*cos(2-2w)^i-f ^wecos(2-2m-c)tfJ, 
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^,=a-im*-. {1 + 3 cos(2-2«i) 6,}, 

we obtain 

p 

—^ = a - Jm*a {1 - 3e cosctf^...- V^« cos (2 - 2m - c) ^J 

/} u 

x{l + 3cos(2-2m)tfJ 
= a + Jm*a {3e + ^me] cosctf,, 

_-, =-f m* {l-4e cosc^j-4wi' cos(2-2m) 0^-^me cos(2-2m-c)tf J 

X sin (2 - 2m) 0^ 
=-fm'{sin(2-2m)^j-2esin(2-2m-c)^j-r^?wesinc^J, 

^ =-a{e 8inc^,+2m* sin (2 -27w) ^j+ ^^Twe sin (2 -2m- c)tfj ; 

therefore 7a-««;73=.**f f^'^ cosctfj, 

the other term is of the fifth order. 

r T 

lyj^ rftf=+f m* cos (2-2m) d^-Zm^e cos [2'2m'C)0^-^rf?e cosctf^, 

(^ ■*'^) ^^ ^^•"" 3m'cos(2-2m) ^j+ ^m'ecos(2-2m-c)tf,; 

^ + w J /^ dO^— ^rri^ae cosc^^, 

and the equation in 2« becomes 

^+w = a{l + (fm'e+ W^'«-|f^'^+V^'e)cosc^,+...l 
= a {1 + (fm'e + ^m'e) cosctf, +...}. 
Assume w = a (1 + 6 cosc^^ +.. .) ; 

therefore ae (1 - c') = (f m'e + ^m'e) a, 

whence c = 1 - f m* - ^m'. 
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104. Hence, to the third order of approximation, 

mean motion of apse _ 1 - c _ f m* + ^^m* 8 + 75m 
mean motion of node ~ ^r — l " ^m^ — /^m' ~ S — Sm ' 

and since m — -^ nearly, we see that the moon's apse 
progredes nearly twice as fast as the node regredes. 

In the case of one of Jupiter's satellites, the periodic 
time round Jupiter is only a few of our days, and the 
periodic time of Jupiter round the sun is 12 of our years, 
therefore w, the ratio of these periods, is very small. 

Hence, the apse of one of Jupiter's satellites progredes 
along Jupiter's ecliptic, with pretty nearly the same velocity 
as the node regredes, assuming these motions to be due to 
the sun's disturbing force; they are, however, principally 
due to the oblateness of the planet. 

Parallactic Inequality. 

105. In carrying on the approximations to a higher order, 
it is found that the expressions for the forces contain terms 
whose argument is the moon's elongation^ i,e. the difference 
of longitude of the sun and moon. These terms appear 
in P and in T of the fourth order^ and since 1 — w, the 
coefficient of d, is near unity, the terms will become of the 
third order, and therefore of considerable importance in the 
values of u and of t We shall work out these expressions 
from the earliest steps ; and we select these terms on account 
of the peculiar use which has been made of them to determine 
the sun's parallax, whence they have received the name of 
parallactic inequality. 

Let us go back to the expressions for the forces 
Arts. (23), (24). In the first place we may remark that 
the terms we are about to compute must be independent 
of 8 ; because any term into the composition of which 8 enters 
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-a 1+ — 7— cos( 



r 
SEG 

r' 



'will necessarily have g in the coefficient of the argnment. 
We may therefore at the outset suppose the orbit to coincide 
with the ecliptic. We have then 

8M* = r" - 2r'MG cos (^ - 6') + MG^, 

(^-5')+^{l+¥co82{^-^)l], 

^cos(^-^') + ^f {|+¥cos2(5-^')}] . 
We must now substitute these values in 

and pick out the terms which have for argument {6 — ff). 
In '«*' (■olus + "ops) t''^'"® ^ill ">* * t^r™ 

3m' J, cos(^-^) = -^ -^-j-^cos(5-^. 

In mVf ojni~ qbs) co8(^ — 5') we shall have the terms 



m 



.{MG'-EG*) 



||co8(^-e') + 1^ cos{0-0')coii2{0-ff)], 



which produce 

m' i pj ■' (1+ V) cos(^-^) = 'i-pr Em'^''^ ^' 

and finally miV (-njjTi - "ops I sin(tf- ^) will produce 



,MG'-EG' 



m 



;V E-M . 



(|_V>)sin(^-^') = |^ J_58ixi(5-5') 
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Therefore the new tenns will he 

Since — (= :j J^^ nearly) is of the second order, these terms 

are of the fourth order. 

106. On the principle of the superposition of small 
disturbances, we may compute the effects of these terms by 
themselves as disturbers of the elliptic path. 

M = a{l + e cos(c5 — a)+ }, 

du • / /I 

-^ = — ae sm(c^ — a) ; 



dd 

therefore 

T du , ^ , , 

TF-8 77: = 0, to the fourth order, 
nu du 

-vps +e*=a, to the second order, 
dGr ' ' 

Substituting in the differential equation for u^ we get 
d 



/: 



^ +w=a 1^1 + - tg5^«-^_^ - cos{(l -7n) e^ ^} + ... I . 

Assume w = a 1+ + ^ ^ w ~ cos{(l-?n)^-/8}+... ; 
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therefore A = j^_ J_^ .. = - if ^ 5 

therefore u=a 1+ -if»i ^ ^ — cos{(l-m) ^— /8}+... . 

107. The corresponding term in the value of will also 
be of the third order, 

h = h[^^ + '^«.|=|^cos{(l-m)^-^}+ ]. 

1 12-8 ^^ IS of the fourth order and will not rise in t ; 
therefore 

■ S=K'^ ■►V^l^f' -1(1 -».)»-« +...]. 

'=1 [«+ + V»f=f J ^{(1 --»)« -« +•••] . 

and ^= 2>^ + — Y^^i — ff"" 8m{(l — ?w)pe — /3}-f...« 

Mayer was the first who applied this term to the deter- 
mination of the sun's parallax, by comparing the analytical 
expression of the coefficient with its value as deduced firom 

observation. The values of m and of -p, and therefore of 

^ — j^ being pretty accurately known, — will be determined, 

that is, the ratio of the sun's parallax to that of the moon ; 
but the moon's parallax is well known ; therefore, also, that 
of the sun can be calculated. The value so obtained by 
Laplace was 8*6"; but Hansen, by the discussion of a great 
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number of modern observations, has found 8'92", a value 
which agrees vpry closely with recent determinations by other 
methods. 

Secular Acceleration. 

108. Halley, about 1693, found, by the comparison of 
ancient and modern eclipses, that the moon's mean revolution 
is now performed in a shorter time than at the epoch of the 
recorded Chaldean and Babylonian eclipses. This pheno*- 
menon, called the secular acceleration of the moorCs mean 
motion^ has not yet been fully accounted for. For a long 
time it was altogether unexplained, but in 1787 Laplace 
gave what seemed to be a satisfactory and complete account 
of it. 

The value of p, Art. (54), on which the length of the 
mean period depends, is found, when the approximation is 
carried to a higher order, to contain the quantity e' the 
excentricity of the earth's orbit. Now, this excentricity is 
undergoing a slow but continual change from the action of 
the planets, and therefore p^ as deduced from observations 
made in different centuries, will have different values. 

The value of p is at present increasing, or the mean 
motion is being accelerated, and it will continue thus to 
increase for a period of immense, but not infinite duration ; 
for, as shewn by Lagrange, the actions of the planets on 
the excentricity of the earth's orbit will be ultimately 
reversed, e' will cease to diminish and begin to increase, and 
consequently p will begin to decrease, and the secular 
acceleration will, so far as this cause is coucemed, become a 
secular retardation. 

It is worthy of remark that the action of the planets on 
the moon, thus transmitted through the earth's orbit, is more 
considerable than their direct action. 
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109. In his investigation, Laplace treated the excentricity 
e' of the earth's orbit as a constant in the diflFerential equa- 
tions, and only considered its variation in the final result. 
This of course greatly simplified the work, but ought to have 
been looked upon only as a first approximation ; yet, strange 
to say, the result so obtained agreed almost perfectly with 
that which had been deduced from a comparison of ancient 
and modem observations. ^^ Malheureusement cette mer* 
veilleuse concordance k laquelle Laplace a At peut-6tre una 
partie de I'^clat de sa plus belle d^couverte, est sensiblement 
alt^r^e par les approximations suivantes."^ 

Some years ago Professor Adams took up the question 
oh initio without any limitation, and introduced the variability 
of e' into the difierential equations themselves. This strictly 
correct process gave a value of the secular variation only 
about half of Laplace's which had agreed so well with 
observation. 

The variation of excentricity does not therefore account for 
the whole observed change, and the mean motion of the moon 
is affected by some other cause or causes at present unknown. 
Perhaps the remaining acceleration is only apparent, and 
arises from a gradual retardation of the earth's diurnal motion. 
The effect would be precisely the same, and such a retar- 
dation may be due to the action of the tides. Whatever be 
the cause, the explanation of this unexplained phenomenon is 
a question worthy of the mathematician's most serious efforts. 

Inequalities depending on the Figure of the Earth. 

110. The earth, not being a perfect sphere, will not at- 
tract as if the whole of its mass were collected at its centre : 
hence, some correction must be introduced to take into 
account this want of sphericity, and some relation must exist 

* Font^ulant, Sygteme du Monde, Supplement, 1860. 
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between tbe oblateness and tbe disturbance it produces. 
Laplace in examining its effect found that it satisfactorily 
explained the introduction of a term in the longitude of the 
moon, which Mayer had discovered by observation, and the 
argument of which is the true longitude of the moon's 
ascending node. 

By a comparison of the observed and theoretical values 
of the coefficient of this term, we may determine the 
oblateness of the earth with as great accuracy as by actual 
measures on the surface. , 

111. By pursuing his investigations, with reference to 
the oblateness, in the expression for the moon's latitude, 
Laplace found that it would there give rise to a term in 
which the argument was the true longitude of the moon. 

This term, which was unsuspected before, will also serve 
to determine the earth's oblateness, and the agreement with 
the result of the preceding is almost perfect, giving the 
compression shsi^ which is about a mean between the 
different values obtained by other methods. 

Perturbations due to Venus. 

112. After the expression for the moon's longitude had 
been obtained by theory, it was found that there was still 
a slight deviation between her calculated and observed places, 
and Biirg, who discovered it by a discussion of the observa- 
tions of Lahire, Flamsteed, Bradley, and Maskelyne, thought 
it could be represented by an inequality whose period would 
be 184 years and coefficient 15"". This was entirely con- 
jectural, and though several attempts were made, it was not 
accounted for by theory. 

♦ Pont^coulant, Systhne du MoruUj vol. iv. 
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About 1848, Professor Hansen, of Seeberg, in Gotha, 
having commenced a revision of the Lunar Theory, found 
two terms, which had hitherto been neglected, due to the 
action of Venus. One of them is direct and arises from 
a ^remarkable numerical relation between the anomalistic 
^ motions of the moon and the sidereal motions of Venus 

* and the earth ; the other is an indirect effect of an inequality 
^of long period in the motions of Venus and the earth, 

* which was discovered some years ago by the Astronomer 

* Royal.'* 

The periods of these two inequalities are extremely long, 
one being 273 and the other 239 years, and their coefficients 
are respectively 27*4" and 23*2". ' These are considerable 
^ quantities in comparison with some of the inequalities already 

* recognised in the moon's motion, and, when applied, they 
*are found to account for the chief, indeed the only 

* remaining, empirical portion of the moon's motion in 

* longitude of any consequence ; so that their discovery may 

* be considered as a practical completion of the Lunar Theory, 
^ at least for* the present astronomical age, and as establishing 
*the entire dominion of the Newtonian Theory and its 

* analytical application over that refractory satellite. 't 

Motion of the Ecliptic, 

113. We have seen. Art. (14), that our plane of reference 
is not a fixed plane, but its change of position is so slow 
that we have been able to neglect it, and it is only when 
the approximation is carried to a higher order that the 
necessity arises for taking account of its motion. 

* Beport to the Annual General Meeting of the Boyal Astionomical Sodetji 
Feb. 11, 1848. 

t Address of Sir John Herschel to the Meeting of the Royal Afltronomical 
Society. 
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It has been found to have an angular velocity, about an 
axis in its own plane, of 48" in a century, and the correction 
thus introduced produces in the latitude of the moon a term 
- CO) cos {0 — 0), 

where a> is the angular velocity of the ecliptic, - the angular 

c 

velocity with which the ascending node of the moon's orbit 

recedes from the instantaneous axis about which the ecliptic 

rotates, ^ the longitude of this axis at time ^, and 6 the 

longitude of the moon at the same instant. 




Let TAm be the position of the ecliptic at time <, 
A the point about which it is turning, TA = ^, 
MN the moon's orbit, JIf the moon, and Mm a perpen- 
dicular to the ecliptic ; Tm = ; Mm = lat. = 13. 
i the inclination of the orbit, and N the longi- 
tude of the node. 
Let APlfm* be the ecliptic after a time hu 
Any point whose longitude is L may be considered as 
moving perpendicularly to the ecliptic with a velocity 
0) sin(i — ^). 

Hence, the point N will move in the direction FN with 
a velocity a> sin (N— 0). And N' will move along PN' with 
a velocity a> sin {N- ^) cott ; 

dN 
therefore -^ = co sin (N^ ^) cot*. 

Again, the point of the ecliptic 90** in advance of N 
will move towards the moon's orbit with a velocity 
a)sin(90° + JV-»; 
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di 
thwefore ^ = — o co8 {N- if>). 

Now, cot i, CO, and — — -, — -^ = — may be considered con- 
stant in integrating ; 
therefore iN^ ecu cos (-W— ^) cot i, 

Si = CO) sin {N" ^), 
and if NM— •^, we have 

^ f ^^ ^^ t XT ^\ 

g^ =. ,=. — , — . cos (N- 6). 

^ cosi sini ^ ^' 

Now, sin 13 = sin i. sin -^ ; 
cos^.S^scost.sin^.Si + sini.cos-^.S^ 

= ccf> {cost sin -^ sin [N- (f>) — cos -^ cos {N— 0)}, 
but cosi.sin'^=cos^sin(^— -W) and cos'^=^co&fi cob[0—N)] 
therefore S0 = — c(o cos (^ — 0). 

The discovery of this term is due to Professor Hansen ; 
its coefficient is extremely small, about V5" ; but, being of 
a totally diflFerent nature from those due to successive 
approximations, it was thought desirable to examine it, and 
the above investigation, which was communicated to me by 
J. C. Adams, Esq.,* will be read with interest on account 
of its elegance. 

With respect to the foregoing investigation, perhaps the following remarks 

dN 
will not be superfluous.: — The value of -,- is not the actual velocity of iV, 

but its velocity relatively to the position of the node as determined when the 

motion of the ecliptic is neglected ; its integral is therefore iN the change of 

longitude due to this motion, and in this integration no constant is added, zero 

dN 
being taken for the mean value. The periodic forms of both -^ and dN shew 

that they oscillate about mean values, the time of a complete oscillation being 
♦ Now Lowndean Professor of Astronomy in the University of Cambridge. 
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that required by em {N-<1>) and cos (N-<t>) to go through their cycle. This 
relative motion of the n,ode is analogous to that of a particle moving in a straight 
line under the action of a force varying directly as the distance. 

Similar remarks apply to the inclination. 

It must also be borne in mind that the result obtained Sfi is not the difference 
between the latitude referred to the actual ecliptic and that referred to a fixed 
plane, but is the difference between the calculated latitude referred to the actual 
echptic on supposition of its being fixed, and the correct latitude referred to the 
same actual ecliptic when its motion is taken into account. 

We may obtain an approximate value of the coefficient 

CO) by substituting for it — , where a is the number of seconds 

through which the ecliptic is deflected in one year = 0'48'S 

and n is the number of years in which the node of the moon's 

2'7r 
orbit makes a complete revolution = 18-6 ; for then, — is 

the angle described by the node in one year; therefore 

wa . ^, ^. n d {N- (l>) • ^ X 

•^-- IS the ratio of a> : — , ' , supposing 9 to remam con- 

stant| which is nearly the case ; therefore 

na 9-3 X 0-48" ^,^„» 

Cft> = — = = 1*42 .* 

27r 314 



* This affords the solution of a problem proposed in the Senate-House 
in the January Examination of 1852. Question 21, Jan, 22. 
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CHAPTER VIII. 

HISTORY OP THE LUNAR PROBLEM BEFORE NEWTON. 

114. The idea which most probably suggested itself to 
the minds of those men who first considered the motion of 
the moon among the stars, was that it described a circle with 
uniform velocity about the earth as a centre ; and this first 
rough result is represented in our equations by neglecting ' 
all small terms and writing =ptj u = a. 

It musty however, have been very soon perceived that 
the actual motion is far from being so simple, and that 
the moon moves with very different velocities at different 
times. 

115. The earliest recorded attempts to take into account 
the irregularites of the moon's motion were made by 
Hipparchus (140 B.C.]. He imagined the moon to move with 
uniform velocity in a circle, of which the earth occupied, 
not the centre, but a point nearer to one side. By a similar 
hypothesis he had accounted for the irregularities in the 
sun's motion, and his success in this led him to apply it also 
to the moon. 

It b clear that, on this supposition, the moon would seem 
to move faster when nearest the earth or in perigee, and 
slower when in apogee, than at any other points of her 
orbit, and thus an apparent unequal motion would be 
produced. 
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EXCENTRIC AND EPICYCLE. 103 

Let BAM be a circle, CA a radius, E a point in AG 
near (7; C5, ED two parallel 
lines making an angle a with 
CA. 

Suppose a body M to describe 
this circle uniformly with an 
angular velocity p, the time being 
reckoned from the instant when 
the body was at B^ and the longi- 
tude as seen from E being reckoned 
from the line ED ; 
therefore DEM = tf, BGM ^pt^ 

AEM = 0''a, ACM^pt-a. 

EG 
Now ym is a small fraction, and if we represent it by e, 

we shall have 

sin if = -p^TTf sin AEM 
CM 

= e sin (^ — a) 

= 6sin(p<+ Jf-a), 

1 —6 COS (pi — a) ' 
this would give Jf, and then ^ by the formula O^pt + M. 

This was called an excentric^ and the value of e was called 
the excentricityj which, for the moon, Hipparchus fixed at 
sin 5^1'. 

116. Another method of considering the motion was by 
means of an epicycle^ which led to the same result 

A small circle P4f, with a radius equal to EG of previous 
figure, has its centre in the circumference of the circle BPD 
(which has the same radius as that of the excentric), and 
moves round E with the uniform angular velocity pj the 
body M being carried in the circumference of the smaller 
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104 LUNAR THEORY. 

circle, the radius PM remaining parallel to itself, or, which 
is the same thing, revolving 
from the radius PE with the 
same angular velocity p^ so 
that the angle EPM equals 
PEA. I 

Now, when the angle AEP 
equals the angle ACM of the 
former figure, it is easily seen 
that the two triangles EPM^ 
ECM are equal, and there- 
fore the distance EM and the angle AEM will be the same 
in both, and the two motions will be identical. 

IIT. The value of e being small, we find, rejecting e*, &c , 
if=6sin (^<-a), 
therefore 6 =pt + e sin {pt - a). 

If we reject terms of the second order in our expression 
for the longitude, and make xj = 1, we get. Art. (55), 

=pt + 2e sin {pt - a),' 
which will be identical with the above if we suppose the 
excentricity of the excentric to be double that of the elliptic 
orbit. 

Ptolemy (a.d. 140) calculated the excentricity of the 
moon's orbit, and found for it the same value as flipparchus, 
viz. 

sin 5** 1' = 1^, nearly. 

The excentricity in the elliptic orbit is, we know, about ^q, 
and i*g and -^ will pretty nearly reconcile the two values of 6 
given above. This shews us, that for a few revolutions the 
moon may be considered as moving in an excentric, and her 
positions in longitude calculated on this supposition will be 
correct to the first order. 



Digitized by 



Google 



EXCENTRIC AND EPICYCLE. 105 

Her distances from the earth will not however agree ; for 
the ratio of the calculated greatest and least distances would 

be - — ^ or J^f , while that of the true ones would be - — ^ 

or 14. 

It would, therefore, have required two different excentrics 
to account for the changes in the moon's longitude and in 
her radius vector. Changes in the latter could not, however, 
be easily observed with the rude instruments the ancients 
possessed, and it was very long before this inconsistency was 
detected. 

118. We have said that the moon's longitude, calculated 
on the hypothesis of an excentric, will be pretty accurate for 
a/ew revolutions. 

The data requisite for this calculation are, the mean 
angular motion of the moon, the position of the apogee, and 
the magnitude of the excentricity. 

But it was known to Hipparchus and to the astronomers of 
his time, that the point of the moon's orbit where she seems 
to move slowest, is constantly changing its position among 
the stars. Now this point is the apogee of Hipparchus's ex- 
centric, and he found that he could very conveniently take 
account of this further change by supposing the excentric 
itself to have an angular motion about the earth in the same 
direction as the moon, so as to make a complete revolution 
in about nine years, or about 3° in each revolution.* 

This motion of the apsidal line follows also from our 
expression for the longitude, as shewn in Art. (75). It is 
there, however, connected with an ellipse instead of an 

♦ On the supposition of an epicycle, this motion of the apse could as easily 
be represented by supposing the radius which connects the moon with the centre 
©f the epicycle to have this uniform angular velocity of about 3° in each revolu- 
tion, and also in the same direction. 
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excentric ; and though the discoveiy that the elliptic is the 
true form of the fundamental orbit was not the next in the 
order of time after those of Hipparchus, yet, as all the irre- 
gularities which were discovered in the intervening seventeen 
centuries are common both to Hipparchus's excentric and to 
Kepler's ellipse, it will be as well for us to consider at once 
this new form of the orbit. 



Elliptic Form of the Orbit 

119. We need not dwell on the steps which led to this 
great and important discovery. Kepler, finding that the 
predicted places of the planet Mars', as given by the circular 
theories then in use, did not always agree with the computed 
ones, sought to reconcile these variances by other combina- 
tions of circular orbits, and after a great number of attempts 
and failures, and eight years of patient investigation, he 
found it necessary to discard the excentrics and epicycles 
altogether, and to adopt some new suppositlpn. An ellipse 
with the sun in the focus was at last his fortunate hypothesis^ 
which was found to give results in accordance with obser- 
vation; and this form of the orbit was, with equal success, 
afterwards extended to the moon: but the departures from 
elliptic motion, due to the disturbing force of the sun, are, 
in the case of the moon, much greater than the disturbances 
of the planet Mars by the other planets. 

In Keplers's hypothesis, then, the earth is to be considered 
as occupying the focus of an ellipse, in the perimeter of which 
the moon is moving, no longer with either uniform linear or 
angular velocity, but in such a manner that the radius vector 
sweeps over equal areas in equal times. 

This agrees with our investigation of the motion of two 
bodies, Art (10). 
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Evection. 



120. The hypothesis of an excentric, whose apse line has 
a progressive motion, as conceived by Hipparchus, served to 
calccdate with considerable accuracy the circumstances oi 
eclipses; and observations of eclipses, requiring no instru- 
ments, were then the only ones which could be made with 
sufficient exactness to test the truth or fallacy of the sup- 
position. 

Ptolemy (a.d. 140) having constructed an instrument, by 
means of which the positions of the moon could be observed 
in other parts of her orbit, found that they sometimes agreed, 
but were more frequently at variance with the calculated 
places ; the greatest amount of error always taking place at 
quadrature and vanishing altogether at syzygy. 

What must, however, have been a source of great per- 
plexity to Ptolemy, when he attempted to investigate the 
law of this new irregularity, was to find that it did not 
return in every quadrature, — in some quadratures it totally 
disappeared, and in others amounted to 2"* 39', which was 
its maximum value. 

By dint of careful comparison of observations, he found 
that the value of this second inequality in quadrature was 
always proportional to that of the first in the same place, 
and was additive or subtractive according as the first was 
so ; and thus, when the first inequality in quadrature was at 
its maximum or 5"* 1', the second increased it to 7** 40', which 
was the case when the apse line happened to be in syzygy 
at the same time.^ 

* It would seem as if Hippaichna had felt the necessity for some fnrther 
modificatiozi of his fiist hjpotheaiSi though he was unable to determine it; for 
there is an observation made by him on the moon in the position here specified 
when the error of his tables would be greatest; and at a time also when she 
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108 LUNAR THEORY. 

But if the apse line was in quadrature at the same time as 
the moon, the second inequality vanished as well as the first. 

The mean value of the two inequalities combined was 
therefore fixed at 6° 2(%'. 

121. To represent this new inequality, which was sub- 
sequently called the Evection^ Ptolemy imagined an excentric 
iji the circumference of which the centre of an epicycle moved 
while the moon moved in the circumference of the epicycle. 

The centre of the excentric and of the epicycle he supposed 
in syzygy at the same time, and both on the same side of 
the earth. 



Cti^- — ^ 

Thus, if E represent the earth, 

8 sun, 

M moon, 

c the centre of the excentric RKTm syzygy, 

5, the centre of the epicycle, would also be in 

^ syzygy. 

Now conceive c, the centre of the excentric, to describe a 

small circle about ^ in a retrograde direction cc'^ while B^ 

the centre of the epicycle, moves in the opposite direction, 



was in the nonagesimal, so that any error of longitude, arising from her yet 
uncertain parallax, would be avoided. Ptolemy, who records the observation, 
employs it to calculate the evection, and obtains a result agreeing with that 
of his own observations. (See Delambre, Ast. Ancienne.) 
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in such a manner that each of the angles 8'Ec\ 8'EE may 
be equal to the synodical motion of the moon, that is, her 
mean angular motion from the sun ; 8ES' being the motion 
'Of the sun in the same time. 

Now we have seen, Art. (116), that the first inequality 
was accounted for by supposing the epicycle RM to move 
into the position nn, r and R being at the same distance 
from Ey and rm parallel to RM^'^ the first inequality being 
the angle rEm, But when the centre of the epicycle is at 
R\ and RM' is parallel to r/w, the inequality becomes R'EM\ 
and we have a second correction or inequality mEM\ 

122. That this hypothesis will account for the phenomena 
observed by Ptolemy, will be readily understood. 

At syzygies, whether conjunction or opposition, the centres 
of the excentric and epicycle are in one line with the earth 
and on the same side of it ; the points r and R coincide, as 
also m and M\ Hence mEM' = 0. 

At quadratures (figs. 1 and 2) c' and R are in a straight 
line on opposite sides of the earth, and therefore R' and r at 

Fig. 1. Fig. 2. 




their furthest distance. If, however, M' and m be at the 
same time in this line, or, in other words, if the apse line 



♦ For simplicity we leave out of consideration the motion of the apse. 
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be in quadratures (6g. 1), the angle mEM' will still be zero, 
or there will be no error in the longitude. But, if the apse 
line is in syzygy (fig. 2), the angle mEM' attains its greatest 
value.* 

Ptolemy, as we have said, found this greatest value to be 
2° 39', the angle mEr being then 5** 1'. 

123. Copernicus (a.d. 1543), having seen that Ptolemy's 
hypothesis gave distances totally at variance with the obser- 
vations on the changes of apparent diameter,t made another 
and a simpler one which accounted equally well for the in- 
equality in longtitude, and was ,at the same time more correct 
in its representation of the distances. 

Let E be the earth, OD an epicycle whose centre G de- 
scribes the circle C CO" about E with the moon's mean 
angular velocity. 




Let CO, a radius of this epicycle, be parallel to the apse 

* If Ptolemy had used the hypothesis of an ezcentric instead of an epicycle 
for the first inequality of the moon, an epicycle would have represented the 
second inequality more simply than his method did. Dr. Whewell's Hittory qfikt 
Jnduaive Sdencesj vol. I., p. 230. 

t See Delambre, Aat. Afodeme, Vol. I., p. 116. Whewell's History qf Inductive 
Sdencest vol. i.» p. 386. 
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line EA^ and about as centre let a second small epicycle 
be described, the radii CO and OMhe so taken that 

CO^OM . _ow .CO+OM . .. 
^ = sm5 l,and jj^ — =sin7 40. 

The radius OM must now be made to revolve from the 
radius OG twice as rapidly as EG moves from E8^ so that 
the 9.ngle GOM may be always double of the angle GE8. 

From this construction, it follows that in syzygies the 
angle GE8 being 0' or 180", the angle GOM is 0' or 360**; 
and therefore G and M are at their nearest distances, as in 
the positions G' and G'" in the figure. Then GM^ GO -. OJf, 
and the angle GEM will range between 0* and 5* 1', the 
greatest value being attained when the apse line is in quad- 
rature. 

When the moon is in quadrature GES^dOf or 270', 
and, therefore, GQM ^ISO"" or 540** and^C and M are at 
their greatest distance apart, as in the position G"] then, 
GM=:i G0+ OMj and the angle GEM will range between 
0° and 7"* 40', the former value when the apse line is itself 
in quadrature, and the latter when it is in syzygy. 

124. Thus the results attained by Ptolemy's construction 
are, as far as the longitudes at syzygies and quadratures are 
concerned, as well represented by that of Copernicus; and 
the variations in the distances of the moon will be far more 
exact, the least apparent diameter being 28' 45" and the 
greatest 37' 33" ; whereas, Ptolemy's would make the greatest 
diameter 1*.* 

The values which modem observations give vary between 
28' 48" and 33' 32". 



♦ Ddambre, Ast. Moderne. 
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125. It will not naw be difficult to shew that the intro»- 
duction of this small epicycle corresponds with that of the 
term ^me sin {(2 — 2wi - c) ;?^ — 2^ + a} in our value of 0. 

For, referring to the preceding figure, we have 

OEM= sin OEM=^smOME 

= ^^in{2. 8EG-AEM) 

= -jY/p sin {2 (moon's mean long. — sun's long.) 

- (moon's triie long. - long, of ap8e)}y- 

and OEM being a small angle whose maximum is 1° 19J', 
we may write moon's mean longitude instead of the true- 
in the argument, and also EG for OE'^ therefore, 

OEM== -^^ sin {2 (moon's mean longitude - sun's longitude) 
— (moon's mean longitude — longitude of apse)}, 

= ^^' sin [2 {pt - {mpt + /3)} - {^« - (l-c)^« + a}] 

= 4770" sin {(2 - 2w - c)pt - 2/3 + a}. 

The value of the coefficient is from modem observations 
found to be 4589-61". 

126. In Art. (77), we have considered the eflfect of this 
second inequality in another light, not simply as a small 
quantity additional to the first or elliptic inequality, but as 
forming a part of this first; and therefore, modifying and 
constantly altering the excentricity and the uniform pro^ 
gression of the apse line. 

BouUiaud (A. D. 1645), by whom the term Evection was 
first applied to the second inequality, seems to hint at some- 
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thing of this kind in the rather obscare explanations of his 
lunar hypothesis, which, never having been accepted, it would 
be useless to give an account of it.* 

In Ptolemy's theory, Art. (121), the evection was the 
result of an apparent increase of the first lunar epicycle 
caused by its approaching the earth at quadratures ; but, in 
this second method, it is the result of an actual change in 
the elements of the elliptic orbit. 

D'Arzachel, an Arabian astronomer, who observed in 
Spain about the year 1080, seems to have discovered the 
unequal motion of the apsides, but his discovery must have 
been lost sight of, for Horrocks, about 1640, re-discovered it 
*in consequence of his attentive observations of the lunar 
^ diameter : he found that when the distance of the sun from 

* the moon's apogee was about 45° or 225*, the apogee was 

* more advanced by 25** than when that distance was about 
*135** or 315^ The apsides, therefore, of the moon's orbit 
^ were sometimes progressive and sometimes regressive, and 

* required an equation of 12** 30', sometimes additive to their 
^ mean place and sometimes subtractive from it.'f 

Horrocks also made the excentricity variable between the 
limits -06686 and -04362. 

The combination of these two suppositions was a means 
of avoiding the introduction of Ptolemy's excentric or the 
second epicycle of Copernicus : their joint effect constitutes 
the evection. 



* Apr^ avoir 6tabli lea monyements et lea 6poques de la lane, Boulliaiid revient 
^ Texplication de r^vection on de la seoonde in6galit6. Si sa throne n'a pas fait 
fortune, le nom da moins est rest^. ^ En m§me temps qne la lone avance sor son 
* odne antonr de la terre, tout le systeme de la lune est deplace ; la terre emportant la 
4ane, lejette loin d'elle rapog6e) et rapproche d'autant le p6rig6e; mais cette 
' ^yection ^ des homes fixees/ Delambre, Hist, de V Ast. Mod.^ torn. ii. p. 157. 

t ^m?^!'^ Astronomical Discoveries of Kepler y p. 307. 
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Variation, 

1 27. After the discovery of the evection by Ptolemy, a 
period of fourteen centuries elapsed before any further ad- 
dition was made to our knowledge of the moon's motions. 
Hipparchus's hypothesis was found sufficient for eclipses, and 
when corrected by Ptolemy's discovery, the agreement be- 
tween the calculated and observed places was found to 
extend also to quadratures; any slight discrepancy being 
attributed to errors of observation or to the imperfection 
of instruments. 

But when Tycho Brah6 (a. d. 1580) with superior instru- 
ments extended the range of his observations to all inter- 
mediate points, he found that another inequality manifested 
itself. Having computed the places of the moon for diflfer- 
ent parts of her orbit and compared them with observation, 
he perceived that she was always in advance of her com- 
puted place from syzygy to quadrature, and behind it from 
quadrature to syzygy ; the maximum of this variation taking 
place in the octants, that is, in the points equally distant 
from syzygy and quadrature. The moon's velocity therefore, 
so far as this inequality was concerned, was greatest at new 
and full moon, and least at the first and third quarters.^ 



♦ *It appears that Mohammed- Aboul-Wefa-al-BoiMdjani, an Axabian astro- 
nomer of the tenth century, who resided at Cairo, and observed at Bagdad in 976, 

* discoyered a third inequality of the moon, in addition to the two expounded by 

* Ptolemy, the equation of the centre and the evection. This third inequality, the 

* variation, is usually supposed to have been discovered by Tycho Brah^, six 

'centures later In an almagest of Aboul-Wefa, a part of which exists in 

' the Royal Library at Pans, after describing the two inequalities of the moon, 
< he has a Section IX., "Of the third anomaly of the moon called Muhazal or 

^ Prosneusis*^ But this discovery of Aboul-Wefa appears to have recited 

^ no notice among his contemporaries and followers ; at least it had been long 
' quite forgotten, when Tycho Brah^ re-<^BCOvered the same lunar inequality.' 
Whewell's Hist, of Inductive ScienceSf vol. I. p. 243. 
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Tycbo fixed the maximum of this inequality at 40' 30". 
The value which results from modem observations is 
39' 30". 

128. Wc have already two epicycles, or one epicycle and 
an excentric, to explain the Grst two inequalities: by the 
introduction of another epicycle or excentric, the variation 
also might have been brought into the system; but Tycho 
adopted a different method:* like Ptolemy, he employed 
an excentric for the evection, but for the first or elliptic 
inequality be employed a couple of epicycles, and this 
complicated combination, which it is needless further to 
describe, represented the change of distance better than 
Ptolemy's. 

To introduce the variation^ he imagined the centre of the 
larger epicycle to librate backwards and forwards on the ex- 
centric, to an extent of 40^' on each side of its mean position ; 
this mean place itself advancing uniformly along the excentric 
with the moon's mean motion in anomaly ; and the libration 
was so adjusted, that the moon was in her mean place at 
syzygy and quadrature, and at her furthest distance from it 
in the octants, the period of a complete. libration being half 
a synodical revolution. 

Annual Equation, 

129. Tycho Brah^ was also the discoverer of the fourth 
inequality, called the annual equation. This was connected 
with the anomalistic motion of the sun, and did not, like the 
previous inequalities, depend on the position of the moon in 
her orbit. 



* For a full description of Tyoho^s hypothesis, see Delambre, Hist, de V Ast, 
Mod. torn. I. p. 162, and An Account of the Astronomical Discoveries of Kepler 
by Robert Small, p. 185. 
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Having calculated the position of the moon corresponding 
to any given time, he found that the observed place was 
behind her computed one while the sun moved from perigee 
to apogee, and before it in the other half year. 

Tycho did not state this distinctly, but he made a cor- 
rection which, though wrong in quantity and applied in an 
indirect manner, shewed that he had seen the necessity and 
understood the law of this inequality. 

He did not try to represent it by any new excentric or 
epicycle, but he increased by (8m. 13s.) sin [sun's anomaly) 
the time which had served to calculate the moon's place;* 
thus assuming that the true place, after that interval, would 
agree with the calculated one. Now, as the moon moves 
through 4' 30" in 8m. 13s., it is clear that adding (8m. 13s.) 
sin (sun'*s anomaly) to the time is the same thing as sub- 
tracting (4' 30") sin [suns anomaly) from the calculated lon- 
gitude, which was therefore the correction virtually intro- 
duced by Tycho.t Modern observations shew the coefficient 
to be 11' 9". 

We have seen. Art. (82), how this inequality may be 
inferred from our equations. 

Bediuitton. 

130. The next inequality in longitude which we have to 
consider is not an inequality in the same sense as the fore- 
going ; that is, it does not arise from any irregularity in the 
motion of the moon herself in her orbit, but simply because 
that orbit is not in the same plane as that in which the 



* That is, the equation of time which he used for the moon dijBEered by that 
quantity from that used for the sun. 

t Horrocks (1639) made the correction in the same manner as Tycho, but so 
increased it that the corresponding coefficient was 11' 51" instead of 4' 30". 
Flamsteed was the first to apply the correction to the longtitude instead of the time. 
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longitudes are reckoned, so that even a regular motion in 
the one would be necessarily irregular when referred to the 
other. Thus if NMn be the moon's orbit and tNm the 
ecliptic, and if M the moon be re- 
ferred to the ecliptic by the great 
circle Mm perpendicular to it, then 
MN and mN are 0", 90°, 180°, 270°, 
and 360° simultaneously, but they 
diflfer for all intermediate values: the diflference between 
them is called the reduction. 

The diflference between the longitude of the node and that 
of the moon in her orbit being known, that is the side NM of 
the right-angled spherical triangle NMm^ and also the angle 
N the inclination of the two orbits, the side Nm may be 
calculated by the rules of spherical trigonometry, and the 
diflference between it and NM^ applied with a proper sign 
to the longitude in the orbit, gives the longitude in the 
ecliptic. 

Tycho was the first to make a table of the reduction 
instead of calculating the spherical triangle. His formula 
was 

reduction = tan'J/ sin2Z/ — \ tan*^/ sin4i/, 

where / is the inclination of the orbit and L the longitude of 
the moon diminished by that of the node. 

The first terra corresponds with the term — JA* sin 2 {gpt — 7) 
of the expression for 6. 

Latitude of the Moon. 

131. That the moon's orbit is inclined to the ecliptic wa& 
known to the earliest astronomers, from the non-recurrence 
of eclipses at every new and full moon; and it was also 
known, sinces the eclipses did not always take place in the 
same parts of the heavens, that the line of nodes represented 
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118 LUNAR THEORY. 

by Nhj in the preceding figure, has a retrograde motion on 
the ecliptic, -W moving towards T. 

Hipparchus fixed the inclination of the moon's orbit to 
the ecliptic at 5', which value he obtained by observing the 
greatest distance at which she passes to the north or south 
of some star known to be in or very near the ecliptic, as for 
instance the bright star Eegulus ; and by comparing the re- 
corded eclipses from the times of the Chaldean astronomers 
down to his own, he found that the line of nodes goes round 
the ecliptic in a retrograde direction in about l8f years. 

This result is indicated in our expression for the value 
of the latitude by the term k sin (g0 — 7), as we have shewn 
Art. (85). 

132. Tycho Brah^ further discovered that the Inclination 
of the lunar orbit to the ecliptic was not a constant quantity 
of S"" as Hipparchus had supposed, but that it had a mean 
value of 5** 8', and ranged through 9' 30" on each side of 
this, the least Inclination 4° .58^' occurring when the node 
was In quadrature, and the greatest 5* 17^' being attained 
when the node was in syzygy.* 

He also found that the retrograde motion of the node was 
not uniform: the mean and the true positions agreed very 

* Mm JouniS) an Arabian astronomer (died a.d. lOOS), whose works were 
translated about 60 years since by Mons. Sedillot, states that the inclination of «the 
moon's orbit had been often observed by Aboul-Kassan-Aly-ben-Amajonr about 
the year 918, and that the results he had obtained were generally greater than 
the 5*^ of Hipparchus, but that they varied considerdbly. 

Ebn Jounis adds, however, that he himself had observed the inclination 
several times and found it 5° 8', which leads us to infer that he always observed 
in similar circumstances, for otherwise a variation of nearly 23' could scarcely 
have escaped him. See Delambre, Hiat. de V Ait. du Moyen Age, p. 1E9. 

The mean value of the inclination is 5° 8' 55*46'',— the extreme valaea are 
40 67' 22" and 6° 20' 6". , 

The mean daily motion of the line of nodes is 3' 10-64:", or one revolution in 
6798*20 days, or 18 y. 218 d. 21 h. 22m. 46 s. 
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CHANGE OF INCLINATION 119 

well when they were in syzygy or quadrature, but they 
were 1° 46' apart In the octants. 

By referring to Art. (87) we shall see that these correc- 
tions, introduced by Tycho Brahe,} correspond to the second 
term of our expression for $. 

Since Hipparchus could observe the moon with accuracy 
only in the eclipses, at which time the node is in or near 
syzygy, we see why he was unable to detect the want of 
uniformity in the motion of the node. 

133i To represent these changes in the position of the 
moon's orbit, Tycho made the following hypothesis. 
" Let ENFhe the ecljptic, -K'its pole, BAG a, small circle, 




having also K for pole and at a distance from it equal 
to 5"* 8'. Then, if we suppose A the pole of the moon's 
orbit to move uniformly in the small circle and in the 
direction BAC^ the node Nj which is at 90' from both A 
and Kj will retrograde uniformly on the ecliptic, and the 
inclination of the two orbits will be constant and equal to AK. 
But instead of supposing the pole of the moon's orbit to 
be at Aj let a small circle abed be described with A as pole 
and a radius of 9' 30'' ; and suppose the pole of the moon's 
orbit to describe this small circle with double the velocity 
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of the node in. its synodical revolution which is accomplished 
in about 346 days, in such a manner that when the node is 
in quadrature the pole may be at a, the nearest point to K^ 
and at c the most distant point when the node comes to 
syzygy, at d in the first and third octants, and at h in the 
second and fourth, so as to describe the small circle in about 
173 days, the centre A of the small circle retrograding mean- . 
while with its unifordi motion. 

By this method of representing the motion, we see that 

when node is in quadrature) . « . f^a=5'8'-9i'=4''58i', 
^ ^tionofthe^^ ^oo^ rA* po,^7* 

syzygy J ^^wtis Kc=^ 8'+9i'=5 n^, 

while at the octants it has its mean value 

Again, with respect to the motion of the node, since N is 
the pole of KaAc^ it follows that when in syzygy and quad- 
rature, the node occupies its mean place; in the first and 
third octants, the pole being at d^ the node is before its 
mean place by the angle dKA=^{^' 30") cosec5" 8' = r 46', 
nearly, and it is as much behind in the second and fourth 
octants. 

So that the whole motion of the node, and the correction 
which Tycho had discovered, were properly represented by this 
hypothesis, which is exactly similar to that which Copernicus 
had imagined to explain the precession of the Equinox. 
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SELECTION OP EXAMINATION QUESTIONS 

FEOM COLLEGE AND SENATE-HOUSE PAPERS AND FROM THE 
MODBRATORSHIP AND FMiLOWSHIP EXAMINATIONS AT 
TRINITY COLLEGE, DUBLIN. 

1. Define the plane of the ecliptic and prove that, as 
seen from the earth, sin A. = ^^^^^ sinX' nearly; X being 
the latitude of the sun and X that of the moon. 

2. Obtain the differential equation of the moon's radius 
vector. 

3. What is the principle on which the successive approxi- 
mations to the moon's path are obtained? 

Compare the coeflScients of the principal term in the cen- 
tral disturbing force with the principal term in the central force. 



4. If in solving the equation 
d'u P T du 



-^ + tt— ,, • 



<%*')i> 



we take as our first approximation M = a {1 + 6 cos(^-a)}, 
we obtain at our higher a{)proximations terms containing 6 in 
their coeflScient. How is this defect avoided by taking 
w = a {1 + e cos (cd — a)} as our first approximation ? 

Upon what principle do we approximate to the value of c ? 

T 
5. ' Shew that in the expansion of ^—s there will be a 

term - ^n^^ sin {(2 - 2w - 2c) 5 - 2/8 + 2a} which would 
rise to the second order in the longitude. 

Shew further that the term with the same argument does 
not appear actwUly in the longitude — assuming 

14 =a [1 + e cos {cQ - a) +. . .+ ^me cos {(2 - 2w - c) 5 - 2/8 + a} 

+...+ Y»«e'cos{(2-2m-2c) d-2)8 + 2a}+...]. 
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122 SELECTION OP EXAMINATION QUESTIONS* 

6. The longitade of the moon contains a term of the form 

A sin{(2-2w-c)p<-2i8-f a} 
where A is of the second order; find how it must have 
arisen and determine A. 

7. Investigate the following expression for the moon's 
longitude as far as the second order, supposing the orbits 
of the moon relative to the earth, and of the earth relative 
to the sun, to be originally circles in the plane of the ecliptic, 

^ =^< + V w* sin {(2 - 27n) pt - 2^}. 

8. Suppose that whatever constants are involved in the 
arguments 5, ^, '^^ &c. in the equation 

i«=-4 + -Bsind+ (7sin0 + i) sin^ + &c. 
are perfectly known, and that observation gives the numerical 
values of a considerable number of longitudes i^, i^, ig, &c., 
corresponding to the known angles d„0j,i/Cj,&c., ^,?<3^»'^? &c., 
^3, <^g, -^3, &c., &c., shew how the constants -4, -d, (J, &c. 
may be numerically determined in the two distinct cases 
where the observations at our disposal are and are not 
unlimited in number. 

9. Explain how this process fails and may be modified 
when two of the terms in question B sin 0, C sin i are nearly 
synchronous in their periods. 

10. Shew that the evection in longitude, viz. 

^^me sin {(2 - 2m - c) ;?« - 2^8 + a}, 
may be represented as the joint eflfect of certain periodic 
changes in the excentricity of the lunar orbit and in the 
mean longitude of its apse. 

11. Assuming the usual notation of the Lunar Theory, 
explain the physical meaning of the following equations : 

5 = * sin(^ — 7), 

s = k s\n{g0^ry)y 

8^k ^m[gd-i)-\-^k sin{(2-2m-5r) ^-2/8 + 7}. 

1 2. Explain the efiect of the term m^a cos {(2-2m) ^-2/8)} 
in the moon's radius vector; find the number of days in 
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tbe period of the resulting Inequality and the ratio of the 
axes of the oval orbit. 

13. Investigate the effect of the annual equation on the 
length of the lunar month. 

Assuming m=-075, <^=-017, sInl3*'30W23345, and the mean 
period of a sidereal revolution of the moon = 27d. 7h. 43m. 12s., 
find the difference between a winter and a summer month. 

14. Shew that the moon's orbit In space Is always concave 
towards the sun. 

Assuming that the earth describes a circle about the sun, 
and the moon a circle about the earth and In the same plane 
with the earth's orbit, compare the curvatures of the moon's 
orbit in space at her perihelion and aphelion. 

15. Prove that If we go to the third order of approxi- 
mation, the motion of the moon's apse In one revolution 
of the moon equals fw' (1 — '^^m) 360." 

P T ,. 

16. The expressions for ^5-5 and tj-s contain terms of 

the form SA cos(^ — d') and -4 sinf^- ff) respectively; com- 
pute the terms resulting from them In the values of u and ff» 

17. Copernicus represented the evection by an epicycle 
superposed on the epicycle which represented the elliptic 
inequality. From the construction given by him shew 
geometrically that the evection vanishes when the moon's 
mean elongation from the sun is half her true anomaly, 
and, In general, varies nearly as the sine of twice the 
difference of these angles. 

18. The secular change of Inclination of the actual ecliptic 
to a fixed ecliptic being a" annually, shew that this will 
give rise to an Inequality In the moon's latitude whose type 

is — jr— a"cos0, nearly; n being the number of years In 

which the moon's nodes make a complete revolution, and 
<l> the difference of longitude of the moon and of the ascending 
node of the actual ecliptic. 
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